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Chapter 1 

Introduction 



This thesis describes how the merging of the two fields of Coulomb drag and meso- 
scopic physics gives rise to interesting physics. The present studies are motivated by 
the fundamental interest in Coulomb coupled systems, but the phenomena that we find 
might as well in the future become important in the context of the increasing packing- 
density in electronic devices on the same chip. Studies of Coulomb drag are also mo- 
tivated in the context of small MOSFETS and drag effects utiUzed for nano-scale 
liquid flow detectors have been suggested recently 

Below we give a brief introduction to the two fields and summarize some of the 
important properties of the electron gas with emphasize on situations where Coulomb 
drag has been studied. Finally, we give an outline for the rest of this thesis. 



1.1 Mesoscopic physics 

Mesoscopic physics has from the field of solid-state physics ^ grown into a rich and 
surprising field itself, see e.g. ^ |[ 0, It deals with systems of often reduced di- 
mensionality at the borderline between the microscopic regime of atoms and molecules 
and the macroscopic world. When the dimensions become shorter than the phase- 
breaking length ^0 (typically on the sub-micron scale at liquid Helium temperature 
for semiconductor heterostructures) the wave nature of the electrons reveals itself fully 
in terms of significant corrections to the classical behavior of the electrons as well as 
completely new phenomena. The laws of classical mechanics do not fully succeed in 
describing the observed phenomena, but on the other hand a full quantum mechanical 
description based on a scaling up of the microscopic description is often impractical 
due to too many degrees of freedom. A successful description is in many cases based 
on concepts from both classical electro-dynamics, quantum mechanics, and statistical 
physics. 

The list of mesoscopic phenomena is long and we shall only mention a few im- 
portant examples; In ballistic wires and quantum point contacts the conductance is 
quantized in units of 2e'^/h |l^] and in disordered wires the sample-to-sample fluc- 
tuations of the conductance have a universal magnitude of e'^/h [O O]. Interestingly, 
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these effects have their counterparts in short (L ^ = hvp/irA) Josephson junctions 
where the wire is contacted by a superconductor at each end. For balHstic junctions the 
critical current is quantized in units of e A / h Jl3[ ] and for disordered wires the sample- 
to-sample fluctuations of the critical current have a universal magnitude of eA/h [Q. 
Another example is the Aharonov-Bohm effect in a small ring connected to leads, see 
e.g. Ref. [|l5|]. Here, a magnetic flux (f> enclosed by the ring causes conductance oscil- 
lations as a function 4> with a period 4>o = h/e{oi higher harmonics 0o / n). Within the 
scattering approach [|l6| the effect can be understood from interference between differ- 
ent Feynman paths circulating the enclosed flux [|l7||. The effect is a clear signature of 
phase-coherence and the damping of the higher harmonics can be attributed to a finite 
phase-coherence length, see e.g. Ref. |]l8[|. 

Not only the electronic degrees of freedom are subject to the bounds of quantum 
mechanics - also the phononic degrees are so. The latter manifests itself in e.g. heat 
conductance quantized in units of T:^k^T/3h in suspended constrictions [|l9[ ^ pT] ]. 
The "universality" is somewhat weaker in the examples with supercurrent and heat 
transport where the "quantum unit" includes either the superconducting pairing poten- 
tial A (material dependent) or the temperature T whereas for the conductance only 
fundamental constants are involved; Planck's constant h and the electron charge e. 

Along with the phase-breaking length i^f, the other important length scale is the 
Thouless length Lt- It may be thought of as the thermal length over which phases 0(e) 
of the quasi-particle states can be considered constant for e in a range of kT around the 
Fermi level. On longer length scales we can not ignore the energy dependence of the 
acquired phase and upon thermal averaging the effect of the strong phase dependence 
may be washed away. 

One of the characteristic things of disordered systems in the mesoscopic regime is 
that sample-to-sample fluctuations are pronounced, though they may not always be uni- 
versal. As the sample dimensions and/or temperature are decreased the phase-breaking 
length and Thouless length will at some point meet the sample dimensions. At this 
point interference effects will modify the classical phenomena and since the interfer- 
ence behavior depends strongly on the particular potential landscape that the electrons 
move in there will be dominating sample-to-sample fluctuations in the presence of even 
weak disorder 



1.2 The electron gas 

Most under graduate text books start with the non-interacting electron gas, see e.g. 
Ref. [^]. Surprisingly, a lot of phenomena in metals and semiconductors can be 
understood from this kind of naive model. The Fermi hquid theory initiated by Landau 
[ ^3[ | explains how a complicated system of interacting fermions may often respond to 
external perturbations just as if it was a non-interacting electron gas. For an overview 
see Ref. [|4| . When studying the excitations of the ground state at least two important 
things emerge; /) the fermionic excitations are of single-particle nature with a spectrum 
that is similar (within e.g. some re-normalization of the mass to an effective mass) to 
that of the non-interacting electron gas and /;) the life time of the excitations diverges 
near the Fermi level. For those reasons these long-lived excitations are often referred 
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to as quasi-particles. 

A large part of the mesoscopic transport phenomena can successfully be accounted 
for by the "scattering approach to quantum transport" p7[ ] — often referred 

to as Landauer-Biittiker formalism — and extensions within the same spirit. For an 
overview of the variety of experiments see Ref. [Q. The formalism applies to the 
non-interacting electron gas so that its usefulness relies on Fermi liquid theory. For a 
review of Landauer-Biittiker formalism see e.g. Refs. [|9[ § ^ 0, 

If all electron systems were Fermi liquids we naively only had to know about the 
non-interacting electron gas and the interaction corrections in terms of the small ratio, 
F ^ 1, between the interaction energy and the kinetic energy [|24|]. However, there 
are various examples of the break-down of Fermi liquid theory or non-Fermi liquid 
behavior, but here we shall only mention a few cases which have also been studied in 
the context of drag. When reducing the dimensionality from three to two (by confining 
the electron gas in one direction) there are already precursors of the break down like 
e.g. logarithmic corrections to the particle-particle scattering rate, — > In T 
[011. 

The break down takes place in the strictly one-dimensional limit. Here the electron 
gas forms a so-called Tomonaga-Luttinger liquid [|3| |3^, 0] — commonly referred 
to as a Luttinger liquid. It is different from the Fermi liquid in the sense that the ex- 
citations are bosonic collective modes rather than quasi-particles. The Luttinger liquid 
has attracted much theoretical interest both because of the analytical progress that can 
be made, but also because of the existence of edge states in the quantum Hall regime 
^ and the discovery of single-walled carbon nanotubes p8[ | which are also be- 
lieved to be one-dimensional conductors Experimentally, it is of course chal- 
lenging to create a one-dimensional system and the most convincing demonstrations 
of Luttinger liquid effects are in the power-law dependence of the density-of-states 
in tunneling studies of quantum-Hall edge states p7[], partly in tunneling studies of 
carbon nanotubes [^, whereas for Raman scattering studies in long semiconductor 
wires some discussion apparently still remain For not too long quantum wires, 
L < hvp/kT, attached to reservoirs through open contacts the Fermi liquid is believed 
to be restored due the strong proximity from the Fermi liquid reservoirs, though some 
debate has existed ^ Fermi liquid theory also breaks down in the extreme 
low-density limit where the screening is reduced so that the energy associated with 
the repulsive Coulomb interactions dominates over the kinetic energy, i.e. F > 1. In 
the ground state the electrons organize themselves into a regular lattice, the so-called 
Wigner lattice, with bosonic excitations which are phonon-like just as for the ordinary 
atomic lattice in solids. The question of a possible metal-insulator transition is as well 
a delicate problem. 

Another interesting state is the quantum-Hall Uquid of the two-dimensional elec- 
trons gas in a strong magnetic field | |4^ , p^ . Depending on the filling of the Lan- 
dau levels the excitations can be so-called composite fermions with fractional 
charge — a quasi-particle with magnetic flux quanta attached at filling factor 
1/ = 1/3,2/5,4/7 etc. Filling factor = 1/2 and other even denominators are special 
cases where there is no quantized Hall effect and where surprisingly the low-energy 
physics (near the Fermi level) is basically that of a Fermi liquid in zero magnetic field 
[B^. The superconducting state of some metals at low temperatures is an example 
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where some other degree of freedom in the system is able to mediate a particle-particle 
interaction which is effectively attractive. When that happens two fermionic particles 
can form a bosonic pair (the so-called Cooper pair) which is allowed to condensate. In 
the conventional low-temperature superconductors — usually referred to as BCS super- 
conductors 1 50 1 — the attractive interaction is mediated by the phonons of the atomic 
lattice, but in the high-Tc cuprates the situation seems to be much more delicate and 
the mechanism is still to be understood. At least for the BCS superconductors the ex- 
citations are known to be of fermionic nature, though they are a kind of superpositions 
of the electron-like and hole-like quasi-particles known from Fermi liquid theory. 

Finally, we will just mention the electron-hole pairs (excitons) formed in semicon- 
ductors consisting of the pairing of an electron excited to the conduction band and the 
oppositely charged hole left behind in the valance band. Such pairs are bosonic and it 
has been speculated that they might undergo condensation. 



1.3 Frictional drag 

When two (or more) electron systems are brought close together many-body interac- 
tions between electrons belonging to each of the two subsystems come into play. One 
of the consequences is the possible momentum transfer between the subsystems even 
in the absence of tunneling mediated charge transfer Driving a current through one 
electron system interaction mediated momentum transfer will show up as an induced 
current in the other layer The induced current is also referred to as the drag current. 
The Coulomb drag effect was first suggested in 1977 by Pogrebinskii and in 1983 
apparently independently by Price [p2[|. However, one could imagine other situations 
such as phonon mediated drag [ p3| , ^4[|55|, drag mediated by van der Waals interac- 
tions [q7|], and plasmon enhanced drag [|58|,p9||. For Coulomb mediated interaction the 
low-temperature behavior is given by a T^-dependence similar to the ordinary rate for 
carrier-carrier scattering, see e.g. Ref. [pO| ], tho ugh corrections emerge when reducing 



the dimensionality from three to two [[31||. Fig. |1.1| shows the temperature dependence 
measured by Gramila et al. [^^. 

The study of Coulomb drag is highly motivated by the fact that the drag conduc- 
tance directly depends on the Coulomb interaction in contrast to e.g the two-probe 
conductance of an open system away from the Coulomb blockade regime, see e.g. Ref. 
iQ. In passing we note that Coulomb coupling also has interest in other contexts, such 
as e.g. capacitive coupling of a mesoscopic conductor to the environmental], 
charge pumping in quantum dots 1 63 1, and also spin polarized transport [ |64| , |65[]. 



Below we will briefly mention a few examples from the development of Coulomb 
drag since the seminal works of Pogrebinskii [ ]5T[ ] and Price [^. For more details on 
Coulomb drag in extended two-layer systems we refer to the review by Rojo [p6|]. 

By now drag has been studied both experimentally and theoretically for quite many 
states of the two electron gases involved. For the bi-layer system the theoretical frame- 
work has been Fermi liquid theory in zero and finite magnetic field |6^ ] though 
still in the integer quantum Hall regime. In the fractional quantum Hall regime the 
framework has been composite fermion theory | |7C| , [7l| ] . 

For quasi one-dimensional drag the studies have been based on Boltzmann equa- 
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Fig. 1.1: Temperature dependence of drag in bi-layer system. The left figure shows 
the approximate T^-dependence and the right figure shows the derivations from this in 
more detail. From Gramila ef aZ. [^^. 



tion approaches ^2, 73 1 within the spirit of Fermi liquid theory and for strictly one- 
dimensional drag between e.g. quantum-Hall edges, semiconductor quantum wires, 
and carbon nanotubes there has been quite some activity based on Luttinger liquid the- 
ory [|74|, [75, 76, 77], ^8|]. A subset of the Luttinger liquid studies has been on drag and 
dissipation-less current in coupled Luttinger rings [ j7^ pO[ ]. 

In passing we mention work on drag between a superconductor and a normal metal 
|]8l|], drag in electron-hole pair systems 1 82 , 83[], and drag between ring-shaped Wigner 
lattices [ |84| ] and between extended metal and Wigner crystals [|5]]. Also the noise 
properties of Coulomb drag has been studied |86|. 

The first experimental studies were on drag between silicon thin-films by Hubner 
and Shockley and drag between a two-dimensional and a three-dimensional elec- 
tron system by Solomon et al. The experimental studies of bi-layer systems were 
initiated by Gramila et al. in the case of drag between two electron gases and by 
Sivan et al. [|8|] for drag between an electron gas and a hole gas. For the subsequent 
studies the main part have been carried out in bi-layer semiconductor systems. As far 
as we know the few exceptions are the very recent studies of coupled quantum wires 

While Coulomb drag of extended bi-layer systems in zero magnetic field is fairly 
well-understood several open questions remain in the context of finite magnetic field as 
well as reduced dimensionality — here we will only mention a few. Recent experiments 
in the fractional quantum Hall regime, filling v — 1/2, have revealed negative drag 
which could be interpreted phenomenologically as drag mediated by a spin-dependent 
interaction [p2|], see also Ref. [|93|]. Also the missing 2kp response [ p4] , p5|] an d the 



low-temperature behavior and deviations from the T -dependence, see Fig. 1.1 



have 



given rise to continued discussion |54 9^ 98, 6^. The same might apply to the 
Luttinger liquid interpretation of recent experiments on quasi one-dimensional wires 
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[ pO| , |9l| ]. Finally, there is the new field of mesoscopic fluctuations and mesoscopic 
Coulomb drag which has only been initiated by a few theoretical papers | ^9| , [TI] , E] 
and which is the focus of this thesis. 

The initial theoretical work was mainly done in the framework of coupled Boltz- 
mann equations [|| H, Hi and first later on a linear-response theory based on 
Kubo formalism was developed by Kamenev and Oreg [ |8l| ] and simultaneously by 
Flensberg et al. [|67[|. Recently drag has also been studied in the framework of non- 
equilibrium Green function formalism [100, 101 1. In the Kubo formalism the DC 



current-current correlation function 11 is calculated by aid of a second order perturba- 
tion expansion in the interaction Ui2- As a result the drag conductivity, 



fT2i(r,r') = ^ (-^) //// dridr2dr'idriC/i2(ri,r2)[/i2W,r^) 

X ^1°° dc.^^^Ai(r',ri,r;,-c.)A2(r,r2,r^,c.), (1.1) 

is expressed in terms of two "triangle" functions A; = —{JiPiPi) and two interaction 
potentials Ui2- This can essentially be thought of as two three-particle Green functions 
connected by two inte ract ion lines which as sugge sted in Ref. [|7| can be illustrated 
by diagram (a) in Fig. L2 The result stated in Eq. ( |l.l| ) will be derived and explained 



in more detail later in this thesis. At this point we n ote that a very similar formalism is 
used in microscopic studies of mechanical friction [ |l02| , |l03| . 

Since there often is some degree of disorder present due to random impurities etc. 
it is useful to know the statistical properties of the drag conductance rather than its spe- 
cific value for a specific disorder configuration. Knowledge about the full distribution 
of the drag conductance is of course the ultimate goal, but often one has to be satisfied 
by the first moments like the mean value and the mean fluctuations. These can be ob- 
tained with the aid of disorder ensemble averaging techniques. For the mean value it 
has become quite common to consider mutually un-correlated (uc) disorder potentials 
of the two subsystems such that for the mean value [p7[ 



(A1A2), =(Ai)(A2). (1.2) 



1.2 



This corresponds to only considering diagrams of the type of diagram (b) in Fig 
where there are no impurity lines crossing between the two triangle functions. This is 
motivated by the situation in bi-layers where the main-source for disorder is the ionized 
dopants located in the outside barriers — the barrier separating the wells is mainly un- 
doped. Disorder scattering in a given well is then mainly due to disorder induced by the 
nearest by outside barrier whereas for the other well the disorder is mainly induc ed by 
the other outside barrier [ |67| ] . However, recently it was suggested by Gornyi etal. [ |l04t | 
to go beyond this assumption and consider the case of mutually correlated (c) disorder 
represented by crossing diagrams like diagram (c). The motivation is that disorder 
induced by the separating barrier will give rise to a common disorder potential for the 
two electron systems. This will especially be the case if the barrier is doped. In reality 
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the disorder is probably only partly correlated. Mutual correlation of the disorder was 



found to increase the mean value compared to mutually un-correlated disorder [104| 




(d) 




(f) 



Fig. 1.2: Diagrammatic representations of current-current correlation functions with 
triangles representing the A,'s and the interaction lines representing the U12S. Dia- 
gram (a) is the bare diagram and diagrams (b) and (c) are the type of diagrams to be 
considered for the mean value. Upon disorder averaging the crossed diagram (c) is 
non-vanishing only for subsystems with mutually correlated disorder Diagrams (d), 
(e), and (f) are the types of diagrams relevant for the fluctuations. Again, upon dis- 
order averaging the crossed diagram (f) is non-vanishing only in the case of mutually 
correlated disorder 

Calculating the fluctuating properties requires the evaluation of two of the bare 
diagrams (a) connected by disorder lines. Studies in this spirit were recently initiated 
by Narozhny and Aleiner who considered mutually un-correlated disorder 



(AiAiAaAa)^^^ = (AiAi)(A2A2), (1.3) 

corresponding to diagrams of type (d). This has subsequently been extended to the 
fractional quantum Hall regime with half-filled {v ~ 1/2) Landau levels by Narozhny 
et al. [^T]]. Finally, in Ref. [E] quasi-ballistic wires were considered where diagrams of 
the type (e) in the presence of mutually un-correlated disorder give the fluctuations to 
leading order in l/kpi with £ being the mean free path. The extension of this limit to 
mutually correlated disorder was studied in Ref. [G]. Including crossed diagrams of the 
same type as (e) it was found that the fluctuations increase by a factor of \/2 compared 
to mutually un-correlated disorder. 
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Fig. 1.3: Diagrammatic 
representation of the aver- 
age of two triangle func- 
tions; (AiAi), (A2A2), 
and (A1A2). 



In general it is a quite difficult job to go beyond the assumption of mutually un- 
corrected disorder for the fluctuations since it requires the average of four triangle 
functions which includes crossing diagrams of the type (f). Furthermore, in the two- 
dimensional case no leading order term can be identified and a systematic infinite sum- 
mation over classes of diagrams has to be performed in order to get a non-divergent 
behavior. Assuming un-correlated disorder one can on the othe r ha nd do with calculat- 
ing the average of two triangle functions like illustrated in Fig. 



1.3 



In fact diagram (c) 

can also be cast into that form and in that sense Gornyi et al. [t04| were very close to 
also calculating fluctuation properties. 

The important thing demonstrated by Narozhny and Aleiner [|9|] for diffusively 
disordered bi-layers was that at temperatures smaller than the Thouless energy, i.e. 
samples smaller than the Thouless length, the sample-to-sample fluctuations become 
pronounced when compared to the mean value. Drag between phase-coherent quasi- 
ballistic one-dimensional wires as well as coupled chaotic quantum dots show the same 
behavior [E]. 

Experimentally, Coulomb drag has not yet been brought into the mesoscopic regime 
with L min(£0, Lt), though experiments on coupled quantum wires are approach- 
ing this limit [ p9| , pO[ |9T| . One of the reasons might be that Coulomb drag increases 
with the sample size and decreases with temperature and since the magnitude of Cou- 
lomb drag is already small in existing samples it is quite a challenge to move into the 
mesoscopic regime. However, the idea of bringing together the fields of Coulomb drag 
and mesoscopic physics is promising and a couple of very recent reports [ ]9^ [7l| , E] 
clearly demonstrate the potential gain of such experimental efforts: Though the meso- 
scopic fluctuations are non-universal they are on the other hand an extreme example 
of mesoscopic fluctuations since they can be of the order of, or even exceed, the mean 
value! 

Two interesting systems are coupled quantum wires and coupled quantum dots. For 
coupled quantum wires there are at least two possibilities based on 2DEGs. They can be 
formed in a bi-layer structure with one wire above the other, each of the wires defined 
by depleting the 2DEGs. The other possibility is to have only a single 2DEG and 
then form two in-plane wires by depleting the 2DEG, see Fig. |L4 



The first situation 
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Fig. 1.4: Different geometries with coupling of two quantum wires. The upper left 
panel shows a scanning electron micro-graph (SEM) image of a GaAs sample with 
wires defined by removing parts of the 2DEG (the dark regions) by a shallow-etch 
technique [p9|]. The upper right panel shows a similar SEM image of a sample where 
the wires are defined by surface Schottky gates [pO|], and the same is the case for the 
sample sketched in the lower right corner [^l]]. The lower left corner shows the simpli- 
fied geometry considered theoretically in Refs. [E, F, G, H]. 

allows for studies of both narrow wires with a few number of modes and also wide 
wires with diffusive transport. The other situations is most suitable for narrow wires 
with a width smaller or comparable to the separation of the wires. For wider wires the 
Coulomb coupling will only be efficient close to the nearby edges of the two wires due 
to screening. 

Another very interesting possibility is to use Nature's own systems such as e.g. 
two nearby carbon nanotubes. In fact, if for multi-walled carbon nanotubes (which 
consist of several concentric carbon nanotubes) one could imagine contacting different 
tubes independently it would in principle allow for a study of drag between concentric 
quantum wires. Whereas this might seem very difficult it is definitely possible that 
only, say, the outermost wall is contacted when performing a two-probe measurement. 
If tunneling between the tubes is absent it means that the current in the other tubes is 
forced to be zero. The electrons in the current-carrying tube will then be subject to a 
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Coulomb friction from the electrons in tubes that do not carry a current and this will 
give rise to additional resistance. Indeed this is very similar to the additional resistance 
due to Coulomb drag between the spin-up and spin-down channels in spin-polarized 
transport [||,|5l|. 




Fig. 1.5: Realization of Coulomb coupled quantum dots which are independently con- 
tacted [105|. The lower right corner shows the generic geometry studied in Refs. [E,F] 
with spatially separated contacts to avoid their otherwise dominating contributions to 
drag. 



The possibility of two Coulomb coupled chaotic quantum dots, see Fig. 1.5, with 
open contacts is an other interesting possibility. For such systems drag will be zero on 
average since when an electron in one dot is scattered by an electron in the other dot it 
can with equal probability leave the dot through either of the two leads. This suggests 
that the fluctuations will be finite. Such systems are of course very difficult to realize, 
but the advances in the fabrication and lithography of semiconductor structures has 
reached a high level of sophistication and very recently Coulomb coupled sub-micron 
quantu m dot s with independent electrical contacts have been fabricated in a two-layer 
system 1 105| even though no drag measurements were performed. 

While the systems with coupled quantum wires [ |8^ , pO| , ^ are examples which are 
close to if not already in the mesoscopic regime (L ^ £^ or L > then the coupled 
quantum dot system [ 105 1 belongs to the mesoscopic regime {\fA < at not too low 



temperatures where there is still sufficiently phase-space for the Coulomb drag to be 
measurable, see Table 1.1. 



1.4 This thesis 



The aim of this work is to formulate a Unear-response theory of Coulomb drag in cou- 
pled mesoscopic systems that otherwise behave as Fermi liquids, i.e. non-interacting 
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Sample 



Dimension Separation Temperature 



Quantum wires [p9|] 


L - 


2 /.im 




^ 200 nm 


< 4.2 K 


Quantum wires [pO|] 


L - 


2 /im 




^ 35 nm 


< 4.2K 


Quantum wires 




2 /.jm 


dr 


--^ 50 nm 


0.2 - IK 


Quantum dots [105| 




- 0.8 /im 


dr 


- 40 mn 





Table 1.1: Properties of the quantum wire samples shown in Fig. 1.4 and the quantum 
dot sample shown in Fig. Note that Coulomb drag measurements have not been 
reported for the quantum dot sample. 



electrons. We closely follow the workof Flensbergef aZ. [p7[|, or alternatively Kamenev 
and Oreg [pT|], on translationally invariant systems. Here, we will however generalize 
the formulation to also include systems with broken translation symmetry and dimen- 
sions smaller than both the phase-breaking length and the Thouless length Lt- 

This is also the approach in the very recent work on diffusive bi-layer systems in 
zero magnetic field and the extension to half-filled Landau levels [ |7l| ] . However, 
here we consider arbitrarily shaped subsystems of non- inte racting electrons with an 
inter-subsystem Coulomb coupling of electrons, see Fig. |1.6| . 

The thesis is organized as follows: Chapters ||, ||, and ^ contain a derivation of 
the formal results for the drag conductance that form the starting point of the cal- 
culations in Refs. [E, F, G, H]. Starting from Kubo formalism the drag conductance 
is calculated with the aid of the Matsubara formalism. In chapter || a discrete for- 
mulation is given which is suitable for numerical implementations. Chapter ^ is on 



V, 




Fig. 1.6: Generic coupled two- 
subsystem geometry. 



V2 
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drag between quantum wires and reports on the numerical and perturbative results in 
Refs. [E,G,H]. Chapter is on drag between chaotic quantum dots and contains the re- 
sults of Refs. [E, F] where the statistical properties are obtain by use of random matrix 
theory. Finally, Chapter gives a brief summary and concluding remarks. 



Chapter 2 

Linear response theory 



2.1 Conductance matrix 

We consider two coupled electron systems in the linear-response limit where the cur- 
rents li are related to the applied voltages Vi through a 2 x 2 conductance matrix, 



h \G2l G22 \V2 



(2.1) 



In the phase-coherent regime one often assumes that the systems consist of non- 
interacting particles and the diagonal parts can then be calculated from the two-probe 
Landauer formula in terms of the scattering matrices. The aim of this work is to study 
the off-diagonal elements in the case of inter-subsystem Coulomb interaction, but with 
electrons of each subsystem being otherwise non-interacting. Furthermore, in order to 
study the pure effect of Coulomb drag we assume that there is no tunneling, i.e. no 
charge transfer between the two electron systems. 

We start from Kubo formalism which can as well be used to derive the Landauer 



formula for the diagonal elements [106, 107, 108|. The Kubo formalism has already 
been used for studying Coulomb drag in extended {L ^ two-layer systems [ pi] , 



and here we shall closely follow the work of Flensberg et al. [^. The main 
difference is that for extended systems the formalism becomes translationally invariant 
(after standard impurity averaging is performed) and thus transformation to the Fourier 
space is convenient. On the other hand this is not possible in the case of mesoscopic 
systems (<C ^0) where a formulation in real space coordinates is required to deal with 
the broken translation symmetry of these systems. 



For simplicity our starting point is a bi-layer system as sketched in Fig. 2.1, but 



the formalism is general valid and more exotic geometries like in Fig. 1.6, are in fact 
included by adding some in-plane confinements. The relative orientation of the two 
subsystems only enters the functional form of the Coulomb interaction between them. 



13 



14 



CHAPTER 2. LINEAR RESPONSE THEORY 



-'1 











' 1 











l2 



Fig. 2.1: Bi-layer system 
with ideal leads. 



2.2 From conductivity to conductance 

In linear response, the induced charge density currents Ji are related to the electric 
fields Ei by the conductivity tensor a through 



Jtir.^t) 'y?fir^,r,,n)E^{r,,t). (2.2) 

Here, Ef{ri, t) = SQ,e''^'' '''~^*\ a and /3 label the spatial directions, and i and j label 
the two subsystems. 

For the study of drag we want to calculate, say, the DC current I2 in subsystem 2 
due to an applied voltage Vi in system 1 so that we can identify G21. From Eq. (2.2) 
we get 




where r2 — (2:2, 2/2) is a position in one of the leads of subsystem 2 with 7/2 being the 
transverse direction. In fact, since we consider the DC limit we can choose the point 
any where due to current conservation, 



V • J, = 0. (2.4) 
The electrical field is given by E'^ (ri) = — V^^$i(ri) with the bias Vi defined by 

Vi = lim $i(ri)- lim $i(ri). (2.5) 

xi — ^ — 00 xi — >oo 

By partial integration we get 



I2 — G22V2 



jj dy2dyi cr^f (r2, ri)Fi 
dy2dri 



^V^^a2t(r-2,ri) 

/3 



*i(r-i), (2.6) 
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where in the second term both xi and X2 are in the leads and we have assumed the two 
leads of subsystem 1 to be identical. In the DC limit the last term is zero due to current 



conservation and comparing to Eq. (2.1) we identify the drag conductance to be given 
by 

G21 =yydy2dyia2T(r2,ri), (2.7) 

where ri (r2) is in one of the leads of subsystem 1 (2). The advantage of this choice is 
that in the leads the asymptotic form of the single-particle scattering states is given in 
terms of plane waves and scattering matrices. In the following we use Kubo formalism 
to establish a quantum mechanical description of 0-21 in mesoscopic structures. 



2.3 Kubo formalism 

In the following we briefly review the Kubo formalism, but in the case of coupling of 
two subsystems. First, we present the Hamiltonian and include the perturbation by an 
external vector potential. Next, we derive the Kubo result for the drag conductivity 

0'21- 

2.3.1 Coupling of two subsystems 

We consider a coupling of the two subsystems described by the Hamiltonian 

= +^ + J^[2, (2.8a) 
where in absence of intra-subsystem interactions we have 

j dr,i,\ri)H,{ri)i,,{r,), (2.8b) 

and for the inter-subsystem interactions 

=^2= / dri / dr2 ■(/'''^(ri)'0'''(r2) C/i2(r-i, r2) ^/'(r2)-0(ri). (2.8c) 



Here, ij] are the field operators. Hi is the Hamiltonian of subsystem i, and Uij is the 
electron-electron interaction between electrons in subsystem i and electrons in subsys- 
tem j. 



2.3.2 Vector potential 



The effect of a vector potential may be included through an additional term. This 
is seen by considering the kinetic part T = ^ {—ih^ + e.Af' of the singe-particle 
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Hamiltonian H = T + V, 



drii;\r)fifi{r)^ /dr-0t(r) 
he 



ilj{r) 



dr ■ 



i2m 



ip\r) V • A(r, t) ip{r) + Vi^(r) A{r, t) ■ Vipir) 



By partial integration we get 

drVi^(r) V • A{r,t)i'{r) 



+ —Air,t)-A{r,t)4'Hr)i'{r)y (2.9) 



drA{r,t) ■ V?/'^(r) ?A(r), (2.10) 



where the surface term at infinity has been omitted. As shown by Sols [107| a more 
careful analysis of the boundary term leads to the same result for the Kubo formula, 
which we will derive in the subsequent section. This means that 



dr?/.^(r)'f ?A(r) = / dr Vi'f (r) TU=o ^(r) 

- y dr (-e)A(r, t) ■ | J(r) + J (r, i)} , (2.1 1) 



where the particle current operator is given by J(r) + J (r, i) with 

J(r) = J!-U^(r)VMr)-\v^Ur)]Mr)] 
2mi I L J J 



J(r,t) 



-A{r,t)ip^r)tp{r). 



(2.12a) 
(2.12b) 

(2.13) 



The full Hamiltonian can therefore be written as Jf + M" where 

^' = - ^ y" dr, (-e)A,(r„ t) • { J,(r,) + J-(r-,, t)} . 

■i 

In the following we "linearize" this term in A and treat it as a perturbation. 

2.3.3 Conductivity 

In the Coulomb gauge where V-A = the electric field is described in terms of 
a time-dependent ve ctor p otential A according \a E ^ —dA/dt corresponding to 
A = —{i/il)E. Eq. (2.13) can thus to first order in E (linear response) be written as 



n 



Y,Jdr,E,ir,,t)-J,{r,). (2.14) 
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In the following we consider the full Hamiltonian + M" where the perturbation 
3^' is switched on at time t = — oo. After the perturbation is switched on, the system 
develops according to the Schrodinger equation 



with 



U{t) 



dt' Jf'{t'^ 



(2.15) 

(2.16a) 
(2.16b) 



where from here on the time-dependent operators (with t as an explicitly stated argu- 
ment) are given in the Heisenberg picture where 



The total charge current can now be written as 



(2.17) 



J,{n,t) = {-e){tp{t)\j,{r,) + j[{r„t)\tlj{t)) (2.18) 
^ (-e)(^(-oo)|C/t(t){j,(r,) + j;(r„t)};7(t)|V(-oo)), 

and in linear response we get 



2 p 

J,(r,,t) {-e){4i~^)\j,ir,,t)\^Pi-^)) + J 

- (VX-cx))|{£;,(r,,0 • / di'e*"(*-*')j,(r„t')}J.(r,;,0|V'(-cx))) 

'J — oo 

+ {ij{-^)\j,{r,,t){E,{r,,t)- I df'e'"(*-*')j,(r„0}|'/'(-«)))}. (2.19) 

J — OO 

Here, we have used that 

U{t) ^ 1 + ^/ dt'g^ydr,£;,(r,,0- J,:(r.,t') 

= 1 + ^ 5] y dr, E,{r,, t) ■ j' dt' e^^^^*"*') J,(r„ t'). (2.20) 
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Let us now consider the current density in, say, the a direction 

- 9 

mil ' ' 



jn{t-t') 



X H)(7^(-cx))|[jr(r„t) , J^{r,,t')]J^{-^))-E^{rj,t), (2.21) 



and comparing to Eq. ([2 .21), we identify afj^ and get 



"^'^ = l^p,ir,)S{n-r,)d,,d^0+'^I\:f^'^in,r„n). (2.22) 



=2 



Here, 



(2.23) 



is the unperturbed particle density and the retarded current-current correlation function 
is given by 

/oo 
d{t-t')e^^^'~''^nf'^ir,,r,,t,t'), (2.24a) 



n: 



ir,,rj,t, t') 



(2.24b) 



where we have used that 



dt'Q{t-t'). 



d{t-t')Q{t-t') .. 



Eq. ( 2.22 ) is the celebrated Kubo formula [109, 35|], but here derived in real space for 
the case of two subsystems so that it in principle applies to all the four elements in 
the conductance matrix i n Eq . (|2.l[ ) . Fo r the diagonal elements the calculation leads 
to the Landauer formula |106, 107, 108|. For the off-diagonal elements we note that 
only the last term in Eq. ( 2.22 ) contributes so that we do not have to worry about the 
usual problem with the divergence of the diamagnetic term in the DC limit, see e.g. 
Ref. JlTol]. 



2.3.4 Finite temperatures 

At finite temperatures, (A) denotes the usual statistical average over the grand-canonical 
ensemble of the quantum mechanical expectation value of the operator A, i.e. 



Tr(V(- 


oo)\e-f^^ A\->P{- 


-oo)) 


Tr (VX- 




oo)) 



(2.25) 
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where /? = 1/kBT (no confusion should be made with the /3 labeling the spatial direc- 
tion!) and the "Kamiltonian" is given by 

= M'- ^N. (2.26) 

For a discussion of the absence of Jf" in the Kamiltonian, see Ref. [^. The trace 
is taken over states \^') with any number (TV) of particles and Q is the grand ther- 
modynamic potential. We shall later make a perturbation expansion in J^u, so that 
these states are eigenstates corresponding to the unperturbed Hamiltonians of the two 
decoupled systems, i.e. \tl') = |V'2)- 



Chapter 3 

Formal calculation of drag 
conductivity 



We are only interested in the transconductivity which is given by the last term in 



Eq. (2.22), 



(f^) = )J^n^f''-(r2,ri,f2). (3.1) 



This is also the starting point of Ref. [g7|. In order to calculate it we apply the Matsub- 
ara formalism, imaginary-time formalism. Turning to the interaction representation we 
make a systematic perturbation expansion in Ui2- After a Fourier transformation we 
utilize the generalized Lehmann representation to obtain an expression for the current- 
current correlation function. 

3.1 Matsubara formalism 

In the Matsubara formalism we write the current-current correlation function as 

n^f(r,r',T,T') = -{i'\f,{j^{r,T)J^{r\T')}\tlj), (3.2) 

where Tr is the usual imaginary-time ordering operator. No confusion with the kinetic 
energy operator T should be made. In the interaction representation we have that 

and due to the usual cancellation of unconnected (unlinked) diagrams [|5]| we get 

n^f(r,r',T,r') = -{i^\fr{mj^{r,T)J^{r',T')}\^), (3.4) 
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where only connected (linked) diagrams are to be considered. Here, = ^ \tp2) 
and 



5(/3) = Tr { cxp 



hp 



In the interaction picture Eq. (2.8c) becomes 



(3.5) 



=^2(t) = J j dridr2 pi{ri,T)Ui2{ri,r2)p2{r2,T), 
where p is the particle density operator. 



(3.6) 



3.2 Perturbation expansion 



Using the expansion of S, Eq. (3.5), we can now calculate IV^i {i' ^r' ,t,t') to any 
order in Ui2- The 0th order contribution Iljf (r, r' , t, r')^°) to the transconductivity 
is obviously vanishing. The 1st order contribution Ilgf (r, r', r, t')*^^^ also vanishes 
in the DC limit [ |67[ | so that the lowest-order contribution is of 2nd order in Uu- We 
will not prove the vanishing of the 1st order term, but just mention that application of 
Fermi's golden rule obviously would give a lowest-order contribution of 2nd order in 
Ui2 too since it involves the absolute square of the matrix element. 



3.2.1 Quadratic contribution 

The 2nd order contribution becomes 



n«f(r,r',r,r') 



A(2) 



1 

2tf 



hp 



dTldT2 



dridr2dridr2 Uuiri, r2)Ui2{r[, r'2) 
X (^|TV{pi(ri,ri)/52(r2,Ti)pi(r'i,r2)/52(ri,r2)J2"(r,r) jf(r',T')}|V). (3.7) 



We can at no cost introduce an extra r-ordering operator TV and since = |?/'i)(X)|i/'2) 
we get 
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Fig. 3.1: Diagram corre- 
sponding to the current- 
current correlation func- 
tion Iljf The hatched 
triangles correspond to the 
functions and A2 . 



1121 ir,r,T,rj ~ 



dri dr2 



dridr2dr'jdr2 



xC/i2(ri,r2)C/i2(r'i,r'2)(^i|f,{pi(ri,ri)pi(r;,T2)jf(r',T')}|^i) 

X (V'2|r4p2(r2,ri)p2(r^,T2)J2"(r,r)}|^2). (3.8) 

Introducing the following three-point correlation function, the so-called "triangle" func- 
tion, 

A?(r, r,r',T',r",T") = -(^,|f.{jf (r, T)A(r', r')/5z(r", r")} |^.), (3.9) 
we get 



n^f(r,r',r,r')(2)=-^ // dridr2 



dridr2dr'jdr2 



X [/i2(ri,r2)C/i2K,r'2)A^*(r',r',ri,Ti,r'i,r2)A^(r,r,r2,ri,r'2,r2), (3.10) 



which is shown diagrammatically in Fig [3.1[ This result can also be found in Ref. 

3.2.2 Three-particle Green function 



From Eq. (3.9) it follows that the triangle function can be expressed in terms of a three- 
particle Green function. We first rewrite the particle current operator 



Ami 



(3.11) 



where {x) = (r, r), and thus 



Af{x,x',x") - ~{fr{jnx)Mx')p,{x")}) 



(3.12) 
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Apart from the order of the field operators this is a three-particle Green function, 

'^z{xi,X2,X3;y3,y2,yi) ^ -{fr{ipi{xi)i)i {x2)^i (xs) (?/3 ) -01 (2/2) V-l ) } ) , 

(3.13) 



so that we get the relation 

fi, . . . , 

(3.14) 



A^(x, x', x") = -^-{V^ - V?}^,;(x, x\ x"- x'\ x', x) 
Zmi 



We note that in order to obtain the correct sequence of field operators, 9 interchanges 
are needed. This gives the additional factor of (—1)^ = —1. 



3.3 Fourier transformation 

— in time, but not in space 

The next step will usually be to introduce a Fourier transformation in both time and 
space as in Ref. However, as already discussed we want to be able to study also 
systems with broken translation symmetry and thus only Fourier transformation in time 
is convenient. The transformation is defined by 



A"fr T r' t' r" r") ^ V 



e 



T_{T — T")—iuJn{r'—r") 



X A^{r,r',r",iLu^,iLu„), (3.15) 

with bosonic Matsubara frequencies, 

Wn = 2mT/h/3, n = 0,±1,±2, . . . . (3.16) 



Applying it to the current-current correlation function, Eq. (3.10), we get 



11/" r'^^ 



dridr2dr[dr2 C/i2(ri, 'r2)C/i2(r'i, 

X A^(r',ri,ri,iw„i,ia;„) A2(r,r2,r2,ia;m/,iw„')|. (3.17) 

Using that J^'^ dr e*('^"^'^">)'^ = h/3Snm it is now easy to carry out two of the Matsub- 
ara sums and get 
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n^f(r,r',T,T')^2) =--^2-^ //// dr,dr■2dr[dr'.,U^2irur2)U^■2{r[,r',) 



2K^ {hpf 

X J2 e"'-*"""'^A{^(r',ri,r'i,icj™,ic^„)A^(r,r2,r^,-iw™,-iu;„). (3.18) 



Finally, we introduce the Fourier transform with respect to the "frequency" of the driv- 
ing field, 

n^f(r,r',zfi„)= / d(r-r')e'^"(---)n^f(r,r',r,r'), (3.19) 
Jo 

with bosonic Matsubara frequencies 

n„ = 2mr/h(3, n = 0, ±1, ±2, . . . . (3.20) 
Applying this transformation we get 

U^^^{r,r',inJ^^ =-^1- ffffdr,dr2dr[dr'2U^2{ri,r2)Uu{r[,r'2) 



2fP h(3 

^ (r', ri,r[, -ir?™, iw„) A2 (r, r2, r'2,i^m, -iw„), (3.21) 



which can rewritten as 

Jl''^P{ry,zn,,f^)^-±-j- nil dr^dr2dr[dr'2U^2{rur2)U,2{r\y2) 



2h'^ h(3 

X ^ A^(r',ri,r'i,-irj„,-icj„) A2(r,r2,r2,irj„,iu;„), (3.22) 



where we have reverted the sum, i.e. iujn — 



3.4 Generalized Lehmann representation 

The aim is now to obtain a formal result within the generalized Lehmann representa- 
tion. First we find the resulting poles and subsequently we carry out the Matsubara 
summation by a contour integral in the complex plane. After obtaining the retarded 
correlation function by analytic continuation we finally take the DC limit. The deriva- 
tion is in principle straight forward but contains many technical steps which may be 



skipped — the final result is given in subsection 3.4.4 

3.4.1 Poles 



We consider the reverse Fourier transform of Eq. (3.15) 
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hp 



Af (r, r', r", ioj^) = / / d(T - r")d(r' - r") 



X e 



SO that the action of the imaginary-time ordering operator in Eq. (3.9) gives 



//N iO„(r-r") 



A^(r, r', r", zf]„, ioj^) = - / d(r - r")e'^ 

Jo 

' d(r'-T")e-"(-'-^")(^.|jf(r,r)p,(r',T')/5.(r",r")|^.) 



+ / d(T'-r")e*'^"(^'-^")(V,|/5,(r',r')Jr(r,T)p,(r",r")|V^) • (3.24) 



We now consider the evaluation within the generaUzed Lehmann representation, see 
e.g. Ref. [Ill], where we denote the eigenstates by | k) with 



(3.25) 



Using Eq. (2.25) this yields 



ft/3 



X / d(r'-r")e-"(^'--")(fc|jr(r,r)p,(r',T')p,(r",T")|A;) 



R/3 



+ / d(T'-r")e-"(--- )(fc|p,(r',T')ir(r,T)p.(r",r")|fc) , (3.26) 



and since the r-dependence of an operator is given by 



(3.27) 



we get 



3.4. GENERALIZED LEHMANN REPRESENTATION 



27 



Afir, r', r", z^i^, zc^„) = -e'^" / d(r - r")e*'^"(^-^ ^ V e"^^-- 



t — t" 



X <^ / d(T'-T")e*""^"-" ^e* 







X (fc|jf(r)e-'*'-^/'"'e'*'-^'/V.(r')e"^'"'/''e^'^''/''/5,(r'')|fc) 
d(r' - r") e('^'=/''+*'^")('''"'^") 



X (fc|p,(r')e-^'"'/''e^'^/''Jf(r)e-^'^/^e^'^"/V.:(r")|fc)|- (3-28) 

Using the closure relation, |m)^m| — 1, twice we now introduce the current and 
particle-density matrix elements and get 

A«(r, r', r", zc.„) = -e'^" / d(r - r")e^""(---") 



krnl 



X (fc|jf(r)|m)(m|p,(r')|0(^|/5»(r")|fc) 



X (fc|p,(r')|TO)(m|jf (r)|?)(?|p,(r")|fc)|. (3.29) 

Performing the inner integrals and using that the Matsubara frequencies are satisfying 
gi/ifi,o„ _ 2 the expression reduces to 



A^r, r\ r", zr!„, za;„) = -e''" / d(r - r") ^ , 



gl)/h+iuir,+in„)(T-T") _ ^({Sk-Sr„)/h+in„){T-T") 



{k\Jt{r)\m){m\p,{r')\l){l\p,{r")\k) 



{Sm — S'l)/h + iuJn 

Sk-<s^)l3^{is^-s,)/h+in„){r-T") _ ^{{g^-s,)/h+iij„+in„){T-T")-^ 
(fc|p.(rO|m)(m|if(r)|;)(;|p,(r'0|fc) | 
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Proceeding in the same way with the remaining integral we get 



Afir,r',r",in„,iLOn) = -e^''Y. 



kml 

f,l3{Sk-gi) _ I 



(A:|jf(r)|m)(m|/5,(r')|0(^|/5.(r-")|fc) 



[S'm - Sl)/h+iu!n 



{Sm - S'l) / h + iVln {Sk- Sl)/h + iuJn+i^r. 

^ {k\p^{r')\m){m\J!^{r)\l){l\j^,{r'')\k) \ 



(3.31) 



This can be further simplified due to the symmetries and interchanging / and k in the 
second term yields 



Ar(r,r',r",zfl„,zo;„) = -e'^^^ 



kral 



{S'k - S'i)/h+ iuJn + ifln {S'k - <S'm)/fl+ irir, 

{k\Jt{r)\m){m\p,ir')\l){l\p,ir")\k) 

{Sm — S'l)/h + iuJn 



{Srn - <S'k)/fl+ i^n {S'l - S'k) / h + ILUn + ifln 

^ {l\Mr')\ni){m\l"{r)\k){k\p,{r'')\l) 

{Si - S'm)/fl + i(^n 



(3.32) 



This result corresponds to Eq. (A2) of Ref. [^. The triangle function is seen to be a 
meromorphic function in the complex z-plane with simple poles at Rc(z) ~ Sm^^k on 
the lines corresponding to Im(z) = — i7„ and Im(z) = and it is analytic elsewhere. 
We also notice that at Re(z) = — Si = Q the pole at iujn = —i^n is removable. 
Since the operators are all Hermitian, the result can be rewritten as 
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i Im(z) 




Fig. 3.2: Contour 
in complex plane 
for the evaluation of 
Matsubara sum. The 
Matsubara frequen- 
cies iuJn are indicated 
by crosses ( X ) and the 
regions where /(z) 
has poles by dashed 
lines. 



X 2 Re' 



E 

kml 



{Sk - + il^n + i^n {Sk - <^m)/fl'+ i^n 

^ (fc|jf(r)|m)(m|p,(r-0|0(/|p,(r'0|fc) 



(3.33) 



from which we see that for = S'l the pole at iw„ = is removable too. We thus 
conclude that Af has no poles coinciding with the Matsubara frequencies icj„ = and 

For later use it is convenient to write A as a function of the three variables ir2„, 
+ iuJn, and iw„, i.e. 



A^r, r', r", zr!„, iw„) A^(r, r', r", zf!,,, zr!„ + iw„, zc^„), (3-34) 
indicating the way in which iri„ and iw„ enter the expression. 



3.4.2 Matsubara summation 



For the evaluation of the Matsubara sum in Eq. (3.22) we consider the function 



/(ifj„,if7„ + z, z) = Af (r',ri,ri, -zf^„, -irj_,z, -z) 



X A2(r,r2,r2,if^„,il7„ + z,z), (3.35) 
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and use the contour 7 shown in Fig. 3.2. Here, S and ^ are infinitesimal small and 



positive and the radius r of the arch is infinite. We consider the integral 

I = f —f{ifln, i^n + Z, z)nB{z) = ^ RCS(/ UB.Zn), (3.36) 

where the Bose function, 

nB{z)^{e^P'-l)-\ (3.37) 
has simple poles at z„ = i{2TT/h0)n ~ iojn so that we get a residue 

Rcs(/n_B,Z„) f{iVtn,iiln + Zn,Zn) Rcs(n_B,Z„) = /(^^",^^n + ^«,^:n) ^ 

hp 

(3.38) 



On the other hand this means that, see Eqs. ( p^ , p5| ). 



1 f (Jiz 

— ^ f{iVLn,i^n + Zn, z„) = j> — /(ifi„, ifin + z, z)nB{z). (3.39) 

The contribution to the contour integral from the arches, 

lim / ^^/(zl7„,za.+r-e*',re'^)«s(r-e^'), (3.40) 

is vanishing since f{i^m + z, z) vanishes in infinity as l/z^. The horizontal parts 
of the contour integral can be expressed as Cauchy principal value integrals so that 

dz f°° dx 

- — : f{iiln,i^n+z, z)nB{z) ^ £P / - — : f{ifln,iil+x+iS, x+iS)nB(x+iS) 
dx 

f{i^n, i^n + X — id, X — i5)nB{x — iS) 

o 27ri 

f°° dx 

+ 0^1 f{iQ.n,x-\-i5,x — i^n + i5)n-Q{x — iVtn + iS) 

J-oc 27rz 

f'°° dx 

+ 0^ I f {iVtn, X — i5 , X — iVln — i5)n^{x — iVln — iS) . (3.41) 

Making the substitution x ^ uj this means that 

1 1 1"^ 

— ^/(il7„,z„) = — ^ / dujUBiuj) 

np ^ ZTTl J-oo 

X |/(ifi„, iQ,n + UJ^UJ + iS) — /{ifln, i^ln + UJ,LU — iS) 

+ fiiVln, LO + 'j(5, LO — iVln) — fiiVln, UJ — iS,LU — if2,i)|, (3.42) 

where we have let (5 ^ in terms with iiln and in the Bose functions n^. Furthermore, 
we have used that nB{x — ifln) = nB{x). 



3.4. GENERALIZED LEHMANN REPRESENTATION 



31 



3.4.3 Analytic continuation 



We now extract the retarded part by analytic continuation [ ]35| , |1 1 1[ ] 

1 1 f°° 

lim V'/(«^^n,^n) = / dwnB(tj) 

in„^n+i(5 2TTI 

X |/(rj + iS, n + LU + i5,LU + iS) - /(fi + i6,Q + ^ + iS, lu - i6) 

+ f{n + i5,uj + iS,u}-n- iS) - f{n + iS,u}-i5,uj-n-iS)y (3.43) 

Shifting the variable lu lu + fl in the 3rd and 4th terms we get 



lim — f(il7„, 



Zn) - / 



iO„ ^12+15 ft/3 

X |nB(cj) [/(ri + + cj + + i5) - + + + - iS)] 
+ 715(0; + f2) [/(fi + O + w + bj-iS)- f{n + i6,n + LU - iS,uj - iS)] |, 



(3.44) 



and rearranging we have that 

X I [71b(lj + ri) — ?ib('^)] + iS,fl + LU + iS,Lu — iS) 
+ n-B{Lu)f{fl + iS,Q + LU + iS,LU + iS) — nB{Lu + Q)f{Q + iS,fl + Lu — iS, lu — iS)^ 



(3.45) 



Using the notation of Ref . [ 67 1 we now introduce the functions 

Af (T, t)A2 (±, ±) = /(±, ±) = f{n + iS,n + LU± iS, LU ± iS), (3.46) 
and get 

iO,.-,f2+i<5 ft/3 ^''^ ' ' 2TTi 

xi^[nBiLu + n)-nB{Lu)] A?(-,+)A^(+,-) 
+ nsMA^^(-, -)A^(+, +) - + +)A^(-, -)}, (3.47) 

corresponding to Eq. (20) in Ref. [|67[|. 
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3.4.4 DC limit 

We now consider the DC limit. 



<^2iir, r') = lim o-jf (r, r', n). 



In Ref. [p7p it was proven that 



lim A^(+, +) ^ lim -) = 0, 



so that Eq. (3.1 ) now gives 



(3.48) 



(3.49) 



X lim , 



dridr2dr'idr2 Ui2{ri,r2)Ui2{r[,r'^) 



A^(+,-), (3.50) 



which simplifies to 



X ^ 



h \ 2h? ) 



duj ■ 



dui 



dridr2dr'idr2 J7i2(ri, r2)C/i2(ri , r'j) 

Af (r', ri, r'l, 0, -uj - id, -uj + i6) 

X A^(r,r2,r'2,0,w + i(5,t^-i(5). (3.51) 



This is the generalization of Eq. (22) in Ref. [ p7p to also include the case of broken 
translation symmetry. 



Chapter 4 

Sub-systems of non-interacting 
particles 



Until this point no assumptions have been made about the unperturbed eigenstates 
of the subsystems; they could be complicated many-particle states or more simple 
states. Assuming that we have Fermi Uquids so that we can consider sub-systems with 
quadratic Kamiltonians 



(4.1) 



we ca n immediately apply Wick's theorem to the three-particle Green function, Eq. 
( 3.14 ). After that we will make an eigenstate expansion and follow the same line as in 
the Lehmann representation. However, the general analysis of poles etc. can be carried 
over to this problem. 



4.1 Wick's theorem 

Applying Wick's theorem [BSD 



where the single-particle Green function is given by 



'^iixuyi) %{xi\y2) ^^(xiiys) 

'^i{x2]yi) '^i{x2;y2) %{x2;y3) 

'^i{x?,]yi) ^j(a;3;y2) %{x3;y3) 



(4.2a) 



(4.2b) 



we get 
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^i{x,x',x";x",x',x) = - '^i{x'\x') 




'Si(x";£) 'Si(x;x") 



'9i(x"\x) 'Si(x'\x") 'Si(x\x') 



»i{x';x) »i{x";x') '»i{x;x") 




'^i{x;x) 



<^i(x';x") 
- '»i{x";x") 

+ <^i{x!';x") 




'^i{x'\i) '^i(x\x') 




%{x'-x:) 




'^i(x;x) 



Fig. 4. 1 : Diagrammatic rep resentation of Wick's theorem applied to the three-body 
Green function, see Eq. (4.3). 



; {^X ^ X {^X ^ X ^ I {^X ^ X )^^2 {^X ^ X )^^2 ( X ^ X ^ 

- ^,(.t; x") + ^.(x; x"). (4.3) 

Only the 2nd and 3rd terms correspond to connected diagrams, see Fig. Ej, so that 



Eq. (3.14) becomes 



(4.4) 



Indeed, since x = x the sequence of the three single-particle propagators are connected 
so that they form a "triangle" as also indicated in Fig. 3.1 . 



4.2. EIGENSTATE EXPANSION 
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4.2 Eigenstate expansion 

In the basis of the eigenstates of the Schrodinger equation, suppressing the subsystem 
index i, 

Jt{r)(bxir) ^ exMr), (4.5) 
the single-particle Green function takes the form jl 1 1[] 



) = E<^A(r-2)</'A(ri)^A(Ti;T2). (4.6) 



A 

Here, the new Green function is given by 



^Al 

X {e(ri - T2) [1 - npiS'x)] - e(T2 - Ti)nFi^x)} , (4.7) 
and np is the Fermi function. 

4.3 Matrix elements 

Introducing the Fourier transform 

^^i^kn) = d (n - T2) e*'="(^i-^^)^^A(Ti; T2) = — -, (4.8) 

Jo ikn-Sx/h 

with fermionic Matsubara frequencies 

ikn = {2n+l)TT/hl3, n = 0, ±1, ±2, . . . , (4.9) 

we get 

^^(a;i;:r2) = ^0Kr2)0A(ri)^^e-*-(^^-^^¥A(*fc„). (4.10) 

A ikn 



Substituting into the first term of Eq. (4.4) we get 



{V^ - V?} ^{x; x')^{x'; x")^ix"; i) 



1 2mi 



X 

Ai A2A3 



X e«-*'^-"»^-^")e('^-"-<)(^'-^")^#Ai(ifc„)^#A,(*A::j^^A3(ifc::), (4.11) 
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with the matrix elements 



P\'x{r) = (l)l,{r)(j)x{r), 



(4.12a) 



and 



(4.12b) 



For the current matrix element to be non-zero there is a special constraint on the two 
states involved in zero magnetic field. Applying the continuity equation V- J(r) = 
to the matrix element Jy\ and replacing the Laplacian by essentially T{r) = J^/f{r) — 
V{r) we get 



= V.JA'A(r-) = (^) |j(<fA' - A)pA'A(r), 



(4.13) 



from which it is seen that the two states in general have to be degenerate, i.e. have the 
same energy. This will be useful later on. 



Proceeding in the same way, the second term of Eq. ( 4.4) becomes 



1 2rai 

x=x " {hj3f ~?r 



Ai A2A3 



X e«-"^")(---")e«-<)(^'-^"'^^Ai(ifc„)^fA.(jfc;)^^A3(«A:;')- (4-14) 



The next step is to Fourier transform with respect to time. 



4.4 Fourier transform 

We consider Eq. ( 3.23| ) and from the first term of Eq. ( |4.4| ) we get 



d(T - r")d(T' - r") e''fi"(^-^")+'''^"(^'-^") 

X {V^ - V^}^(x;x')^(a;';a;")^(x";i) 

X — X 

= ris -T-y^^'^'lr, i"',r",zA:„,iA:„ + iujn.ikn ~ i^^n), (4.15) 
tip h ^ 



where we have introduced the function 
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Ai A2 A3 

X ^'^x^{ik^WxA^kn + iuJnWxSkn " i^n)- (4.16) 
The second term is treated similarly and putting things together we arrive at 

A"(r, r', r", zfi„, zri„ + iuj^, iut^) = 

hp 

X ^ |A'(r,r',r",ifc„,ifc„ + iu!n,ik„ - i^n) 

+ K{r,r",r',ikn,ikn - iw„ - iVln^ikn - (4.17) 

4.5 Matsubara summation 

For the evaluation of the Matsubara sum in the function K we use the identity 

A-^B-^ = {B - A)-^{A-^ - B-^) (4.18) 

so that 



fXi {iknWx2 {ikn + il^nW\3 [ikn - i^n) 



X <^xAikn + ioJn) ['AAikn) - ^xMn ~ i^n)] ■ (4.19) 



Repeating this "trick" we get 



11 -\- iuJiij^x.^(^ikji i^^ji^ 



hp 

_ 1 [ 7^ Eifc„ [^A2 [ikn + iuJn) - ^Ai (»fc»)] 



(4.20) 



(A2 - (^Aa)/^ - i^n - i^n 

To evaluate (ikn) we consider the contour 7 in Fig. 4.2 and the integral 



- f ^Jiz)nFiz) = - ^ Res(/ "f, = ^ E /(^")' ^^.21) 
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where f{z) = (z) and 



nj.(z) = (e''^^ + l) \ {A.22) 



, -1 



is the Fermi function. 

Since the contribution from the arches is zero and the two horizontal parts cancel 
only the part encircling the pole z = Sx^/hof f{z) contributes. Therefore 

|^/(^)"^^(^) = -"^'(^) ^) = -M^)^ (4-23) 

so that, removing the h in the definition of ri^, 

^ ^^Ai(«fcn) = nFi<^x,). (4.24) 

Shifting the pole we in the same way find that 

(4.25) 

In this way we now get 



hp 



(4.26) 



4.6. ANALYTIC CONTINUATION 
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where we have introduced the susceptibihty 



XAA' 



lLUn + iS'x-<^X')/h' 



(4.27) 



Returning to Eq. (4.16) this means that 



ik„ 



J^3M{r)px^M{r')px2M{r") 



XXi ,X2 (il^n) - XX3,X2 ji^n + i^n) 

{Sx,-s^,)lh-in^ 



(4.28) 



and treating the other term similarly Eq. (4.17) becomes 



A"(r,r',r",zf2„,zr2„ + iw„,za;„) = - ^ 



1 



Al A2 A3 



X { Jx3Xi{r)pxiX2{r')px2X3{r")[xxiX2{'i^^n) - XXiX2i.it^n + i^n)] 



[Jx3Xi{r)px,X2ir')pX2X3{r")]* [xAi A2 (-i^^n) - XAiAaC^i^^n - i^n)] >■ (4.29) 



Here, we have used that JA3A1 (r) is non-zero for S'x^ ~ only, see Eq. (1-.13) 



4.6 Analytic continuation 

At this stage we should in principle carry out a standard analytic continuation. How- 
ever, this has already been done within the generalized Lehmann representation and 
we can immediately take over the results obtained for A( — , +) and A(H — ), see Eqs. 
(|3.46|,^3o|,|3.5l|). The first one 



A(-, +) = A(r, r', r", iS, -u - iS, -cu + iS), 



becomes 



A(-,+) = - E 



Ai A2 A3 



{<^x^-<^X3)/fi-^S 



X \ JX3X1 ir)px,X2 (.r')px2X3 (r") [xaiAs {-^ + iS) - XX3X2 - i^)] 



- [Jx3X^ {r)px^X2 {r')px2X3 [r")]* [XX3X2 - iS) - XAiA^ + i5)\ \ , (4.30) 
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or in terms of the retarded susceptibility 

"-fd('^a) - "fd(<^a') 



Xxx' (t^) 



(4.31) 



we get 



A(-,+) = -27r;i WA3Ai(r)pAiA,(r')pA,A3(r-")Im[x^,AiH] 
A1A2A3 I 

- [JA3A,(r)pA,A.(r')pA.A3(r")]*Ini [x^a,a.(^)] - '^Xs)- (432) 

Here, we have used that (?Ai = <^A3 to replace [((#Ai — S'\^)/h — iS\~^ by iTrh6{S'x-^ — 



S'xs), see Eq. (4.13). Since 



Im[xL'(^)] = -TTh[nF{<^x) - npi^x')] ^Sx^ A'), (4-33) 
we finally get 



A(-,+) = 2^2^2 ^ |jA3A,(r)pAiA.(r')pA.A3(r")'5(^A. -^Ai +M 

Ai A2A3 

+ \Jx,xAr)px,xAr')px,xy)X KSx, - Sx, - hw)] 

X [^f (^A.) - npiS).,)] <5(<?Ai - <Cx3)- (4.34) 

In the same way we find that 

A(+,-) = A(r,r',r",^,5,w + ^(5,u;-^(5) = -A(-,+) . (4.35) 

UJ — y — uj 

4.7 DC drag conductance 

In the following we give the final result for the DC drag conductance, in terms of both 
the matrix elements introduced in Sec. 4.3 and also in terms of spectral functions. 



4.7.1 Matrix element formulation 



From Eqs. (2.7, 3.51 ) it now follows that 



G21 



X ^ 



^^^A,(.,.-,<)A.(-.,.-,.-)^ (4.36) 



2fcTsinli (Sw/2fcr) 
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where we have defined Ai{uj, r^, r") = ~ih ^ J dyi A(— , +) which is given by 

A1A2A3 

X 5(^2 - + M + [-^A3AiPAiA2('^DPA2A3('''/)]*'5('C\2 - ^1 - | 

X - ni.(AJ] ^(Ai - ^3)- (4.37) 

Here, /A3A1 = / dj/i JA3A1 (t";) has been introduced. There is no spatial dependence of 
the current matrix element /A3A1 due two current conservation. To arrive at this result 
we have used that 



1\ dnB{i^) 
2 



duo 8fcTsinh^(?kj/2fcr)' 
and assumed spin-degeneracy so that 

A1A2A3 A1A2A3 



(4.38) 



(4.39) 



Eqs. ( [4.3q , |4.37[ ) is the generalization of the existing theory of Coulomb drag to also 
include the mesoscopic regime with broken translation symmetry. This formulation 
constitutes t he fo rmal starting point in Refs. [E, F, G, H]. We will end this section by 
writing Eq. (4.37) in terms of spectral functions. 



4.7.2 Spectral function formulation 

The aim is to write the triangle function in terms of the spectral function 



^£(^1,^2) 



27r^0^(ri)</.A(r2)J(e-eA). (4.40) 



To do that we first write out the matrix elements and re-arrange the order of the eigen- 
functions 



JA3Ai(r)pAiA2(T'')PA2A3(^") 



2mi J r^r ^ 



X (r')0A, (rWx, (r")0A2 {r' Wx, ir)<Pxs (r"), (4.41) 



and in Eq. (4-. 37) we thus now have three pairs of eigenfunctions, three sums, but only 
two delta functions. We can at no cost get the third delta function by introducing 
1 = / d(f (5((f - S'x2) into Eq. (0%. Putting things together this gives 
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X {£/g+n^{r'^,r,)j!/g{r'^,r'^)£/g+n^ir„r'^') [npiS') - npiS" + fuj)] 
- [■i^s-Ur-,r,)j^dr-,r'i).s^^^nUn.r'^)]* K(^) - npir^ ~ hw)] }. (4.42) 

In the first term we make the shift S" ^ S" — Huj and in the second term + hu 

and taking the Hmit f ; r-i we get 

(4.43) 

This resuh is the starting point of Ref. [G]. It turns out to extremely useful for numerical 
implementations since efficient methods for obtaining the retarded Green function 
in mesoscopic structures has already been developed. The original purpose of these 
methods was to obtain the diagonal part of the conductance matrix via the Fisher-Lee 



relation [106| and other properties like the density-of-states which can be formulated 



in terms of scattering matrices |112]. For the diagonal conductance it is sufficient to 
know the amplitude of propagating from one boundary to the other (one lead to the 
other) whereas for the drag conductance the full Green function is needed since the 
drag is mediated by Coulomb interactions between particles located anywhere in the 
two subsystems. 



4.8 Low temperature expansion 

At low temperatures fceT <C /i we expect the same temperature dependence of the drag 
conductance as for the ordinary electron-electron scattering rate since the mechanism 
is the same. The available phase space for electron-electron scattering of course goes 
to zero at very low temperatures and thus the Coulomb drag conductance decreases 
with decreasing temperature and vanishes at T = 0. At low temperatures the two 
Pauli blocking factors entering the electron-electron scattering rate gives rise to a 
dependence . As as we shall see we also get at quadratic temperature dependence 
of the drag conductance. 



We start by noting that the sinh factor in Eq. (4.36) cuts of the frequency inte 



gration and at low temperatures we can expand the triangle function to lowest order in 



uj. This gives A oc w and since in Eq. (4.36) 



nd. . ^(fcT)^ (4.44) 



4.8. LOW TEMPERATURE EXPANSION 



43 



we get G21 oc r^. In the following we make the systematic expansion of the triangle 
function that gives rise to this result. 

4.8.1 Matrix element formulation 

First we note that 



X lim 



dS' 



+ 0(^2) = hjS{S" - S')S{^) + 0{uj^), (4.45) 



and substituting into Eq. (4.37) we get 



Ai A2 A3 



As expected the low temperature behavior is governed by the states at the Fermi level 
where all the matrix elements are evaluated since 6{S') = S{£ — ep)- This expression 
is the starting point for the studies of chaotic quantum dots in Refs. [E, F] and the 
perturbative studies of weakly disordered quantum wires in Refs. [E, G, H]. 

4.8.2 Spectral function formulation 

Similarly we note that 



huj^/g lim 



OS" 



+ 0{uj^) = hLU.S!/sSi(?) + 0{uj^), (4.47) 



and substituting into Eq. (4.43) and performing the energy integral we now get 



Zi: Zm L 

X j dy,[^,,(r:.,n)a,Xep(n,0-^s.(n,rn5,,^e.K,r,)]}, (4.48) 

with all spectral functions evaluated at the Fermi level. This expression is the starting 
point for the studies of quasi-ballistic ID wires in Refs. [E, G, H] where it is studied 
numerically by mapping it to a lattice. 
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4.8.3 Scattering basis formulation 

We consider Eq. ( [4.46| ) with the matrix elements calculated in the basis of scattering 
states, 

>r,^/'2",---'0^'---), (4-49) 

where tpl^ is a scattering state incident in sub-band (mode) n from either the left lead 
(I = +) or right lead {I = — ), and is the number of propagating modes. 

The current is conserved and thus we are free to evaluate the current matrix at any 
point we like. In particular we choose to do it within the leads where the asymptotic 
behavior of the scattering states are given in terms of plane waves and the scattering 
amplitudes entering the 2N x 2N unitary scattering matrix S, which has the block 
form 



f,), (4.50) 



with the blocks being N x N matrices. The states incident from the left lead have the 
form 

f 4't{r) +J2n'rn'n(l)n>{r) , X < 0, 

^+(r) = <^ (4.51a) 
and similarly for the states incident from the right 



^-(r) = <; (4.51b) 



Here 



,±/ ^ _ exp(±iA:„a;)T„(y) 



^^{r) - -^v^-p'-'^^^ , (4.52) 



where the transverse wave functions are orthonormal 



j AyT*^{v)Tr.{y)^5^n. (4.53) 



The normalization in Eq. (4.52) makes all modes carry the same current. 

The previously introduced quantum number is expressed as A = (?/, where 
T] = {I, n). Since 

^ = h^kl/2m + en- ^Ji, (4.54) 
with e„ being the sub-band energy, the summations would normally be expressed as 
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n 



However, the fc„ ^ = [\/k^^/Ki\ ^ factor in [dS/dkn) = m /h^kn is already in- 
cluded in the normalization of the wave functions, see Eq. ( p2| ), and we get 



^ — > ^ df^ Q{(ff + fi ~ En) [ scattering states! ]. (4.55) 



From Eq. (4.46) it now follows that, suppressing the subsystem index i. 



Aiu;,r',r") = {2nfh(^^y /jt^ImTr {/p(r')p(r")}, (4.56) 

with all matrices evaluated at the Fermi level. Since the matrices are Hermitian we can 
now also write the result as 



A(a;,r',r") = l(2^)2r.(^)' hu Tr {l[pir'); p{r")] _} , (4.57) 

where the commutator of two matrices M and N is defined in the usual way [M ■,N]- = 
MM - NM. 

Current matrix 

— cancellation of velocity and density-of-states 

We now calculate the current matrix 



/ = 



I++ I+- 
/-+ I— 



(4.58) 



Here, it is most convenient, but not necessary, to consider the right lead (x > L) for 

/++ and the left lead {x < 0) for / This gives I^^ — {h/m)tH and similarly 

= —{h/m)t'^t' . Though slightly more complicated the off-diagonal matrices 
/ ^ = /|_ follow in the same way so that 

h ( tH -rH'\ 

Note that the dimension of / changes to that of h/m with our deffiniton of scatter- 
ing states. By a straight forward matrix multiplication and use of the unitarity of the 
scattering matrix it is easy to show that the matrix can be written in the compact form 
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I=^[t'-S^t^S], (4.60) 



where t is a block matrix with the form of the 3rd Pauh matrix 



-1 



(4.61) 



where the blocks are x iV matrices. As seen the velocity hkn/m is canceled by the 
density-of-states, i.e. [\/k^\/k^^^ from the normalization of the wave function, see 
Eq. ( 4.52 ). This is a central point in the derivation of the Landauer-Biittiker formalism 
for the conductance - essentially the diagonal parts of / - but of course also in cal- 
culation of the off-diagonal elements which are important in the context of Coulomb 
drag. 



The compact notation in Eq. (4-. 60) will be useful when studying drag between 
chaotic quantum dots since in that case standard random matrix theory applies 
immediately to the scattering matrix S and thereby to /. 



Particle-density matrix 
— van Hove singularities 

For the particle-density matrix p it is obvious that there is no cancellation of the density- 
of-states, i.e. 



{Pn'(r)}„„, ={^i(r)}*Vi(r), 



(4.62) 



will have a singular behavior — van Hove singularities — if the Fermi level is at 
the onset of mode n or n' due to the normalizations l/VX^ and Xjyfk^i of the wave 
functions, see Eq. (4.52). In fact this will be the case for any mode for an arbitrary 



scattering matrix that allows for mode-mixing, see the the sum over n' in Eqs. (4.51 ) 



The apparent divergence of the drag at the onset of new modes is of course a failure of 
the low-temperature expansio n and the built-in thermal averaging in a small window 
around the Fermi level in Eq. ( 4.37 ) removes this artifact. Though non-divergent there 
however will still be a peak in the drag conductance at the onset of new modes. This 
was also found in a theoretical study of drag between ballistic wires with a multiple 
number of sub-bands [O], but here it is proven for an arbitrary scattering potential. 



Capacitive and point-like coupling 
— no contribution to Coulomb drag 



Some very general consequences follow directly from Eq. ( 4.57). First we notice that 



the general property A((jj, r', r') = follows directly from the commutator. However, 
even more importantly it is seen that if we in Eq. ( 4.36 ) consider a capacitive (constant) 
coupling. 



U\2{i'\,r'^ ~ Ui2 = constant, 



(4.63) 
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then G21 = since it is easily realized that 



J I dr'dr"A(tj,r',r") cx fiwTi {l[ j dr'p(r') ; j Ar"p{r")\_] = 0. (4.64) 

This conclusion also follows in concrete calculations on simple systems like e.g. bal- 
listic wires, but for an arbitrary disorder potential this property is difficult to show 



explicitly without Eq. ( 4.57). The physical reason is that the electrons do not repel 



(or attract) each other, i.e. when two electrons interact they might stay very close to 
each other as well as far apart. Both situations are equally favorable with respect to the 
interaction energy. 

If the two subsystems interact in a single point R, i.e. 



U,2{r[,r'^)^5{r[~R)S{r'^~R), 



(4.65) 



it again follows directly from Eqs. ( 4.36 , 4.57 ) that G21 = 0. We note that this does 
not exclude neither a short ranged interaction like 



Uu{r[,r'^)^5{r[-r'^), 



(4.66) 



from contributing to Coulomb drag nor a set of point-like couplings. 



Chapter 5 

The discrete problem 



The numerical method that we will develop in the following is — as is usual the case 
— based on a mapping of the continuous problem onto a discrete lattice problem. For 
each subsystem we choose a finite set of spatial points where the needed functions and 
quantities are calculated and evaluated. In general the set of points may form a regular, 
say, square lattice for two-dimensional problems or it might have some irregular and 
even three-dimensional form. In this work the lattice only has a mathematical meaning 
and the lattice constant is to be chosen small enough that the convergence towards the 
continuous problem is ensured. However, the lattice could also be the underlying phys- 
ical lattice of the atoms like it would be the case if we were going to make a detailed 
study of drag between molecular wires like e.g. carbon nanotubes. In what follows we 



will consider two square lattices as illustrated in Fig. 5.1 which is the discrete version 



of Fig. 2.1 The way we number the lattice points is not crucial and here we will for 



simplicity choose the one indicated in Fig. 5.1 though at the end the overall formulation 



will not depend on the particular way of numbering. 

5.1 Trace formula for drag conductance 



We consider Eq. ( }4.36[ ), where both the triangle function A and the interaction Uu 
are functions of two spatial variables and can thus be represented as two-dimensional 
matrices with elements 



{Ul2}mi7n2 = C^12(^•l("^-l),f2("^2)), (5.1) 

{A,([j)}„^„^ = A,(w,r,(m,;),ri(m-)). (5.2) 
Approximating the spatial integrals by sums over the lattice sites, 

/ dr ^ (V/M) ^, ViD = L, V2D = WL, 
49 
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m=l 




m=M=MyMx 2Mx^^ 




a 



L=My-a 




Fig. 5.1: Bi-layer system of width W and length L with each 2DEG represented by an 
Nx X Ny square lattice with lattice constant a. 



we get 



G21 = ^ X ^ 
h 



Tr{[/2iAi(c^)L/i2A2(a;)} 
2kTsmh^{njjj/2kT) ' 



(5.3) 



where we have used that Ai( 



Ai{uj, r", r'J and introduced U21 = Ui 



The "missing" factor (V/M)^ will be incorporated into the spectral function so that 
will have the dimension of energy. 

For simplicity we will assume the same number of sites in the two subsystems 
so that all matrices are M x M matrices instead of the general case with A/s being 
Ali X Mi matrices and U12 being an Mi x A/2 matrix. 

The one-dimensional case {My = 1) is most illustrative and the generalization to a 
finite number of modes will be given afterwards. 



5.2 One-dimensional formulation 

For the one-dimensional problem we label the points by in from 1 to M and formulate 
the problem in terms of M x M matrices. 



5.2.1 Coulomb coupling 

For the Coulomb coupling the interaction between site m in one wire and site m' in the 
other wire is given by 

{?/i2},„,„, = Ui2{x[ ma,x'2 m'a). (5.4) 

In most cases C/12 will be symmetric, but e.g. a spatially dependent screening would 
break this symmetry. 
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5.2.2 Triangle function 

For the triangle function we similarly have that 



{^'(^)}mm' = ^^i^^^'i rna,x'^ m'a). (5.5) 



It turns out to be most convenient to work with spectral functions, Eq. (4.43), since 
the spectral function ^ is easy to obtain numerically whereas the eigenfunctions only 
follow implicitly. 

Similarly to the transformation of integrals into sums we now need to transform the 
spatial derivative into finite differences, see e.g. Refs. For some function F this 

means that 



OxFlx = ma) — > — — ^ (5.6) 

a 

which of course only becomes a good approximation when the lattice constant a be- 
comes sufficiently small. We will return to a more detailed criterion for that later. 

Due to current conservation the point at which the derivative is taken can be any- 
where. Applying Eq. (5^) to Eq. ( 4.43 ) at a; = ma we get 

L L ^ J mm I- J m,m-f-l L 1 vam ^ J m+l,m J 

X [npis — hcu) - nF{e)] - ( a; ^ -lo )*. (5.7) 

In the square brackets two out of four terms have canceled when introducing the finite 
differences approximation. Here, the spectral function matrix has matrix elements 

{■^e}j^j-^, = a X £/^{x ma, x' m'a). (5.8) 

The normalization by the lattice constant a gives the dimension of energy — it is 
a matter of choice which turns out to be convenient. 

Due to current conservation we are free to sum over rfi and divide by the number 
of lattice points corresponding to integrating along the wire and dividing by the length. 
This sum can be incorporated by a matrix product. Introducing the matrix A with 
elements 

-L r r r f ±Om.,m'±l 



M _ 

rh—1 



^ / ^ ^ni , m + 1 f^m 771 ' rh ^rii + 1 , rn ' — -j^ ; (5 ■ 9 ) 



we get 



X [npie — hhj) — nF{e)\ — [ uj ^ -uj )*, (5.10) 
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where (g) denotes an element-by-element multiplication, {X ® Y}nm = XnmY7im- In 



Eq. (5.9) we only sum over the first M — 1 sites since Eq. (p.q) can not be applied to 



the site AI if we insist on having M x M matrices. 

5.2.3 Spectral function 

The spectral function £/ is obtained from the retarded Green function since from 



Eq. (4.40) 



i4 = i - = i - . (5.11) 

The problem is thus reduced to that of finding a method which provides the retarded 
Green function for any given scattering potential. 

5.2.4 Hamiltonian 

To get the Green function we first map the continuous Hamiltonian 

.ye ^-^&l + v{x), (5.12) 

onto a lattice. We do that by applying the finite differences approximation to the Lapla- 



cian. Eq. (5.6) is not a unique definition and we might as well have considered the other 
neighboring point m — 1. Using both neighboring points for the Laplacian we get a 
symmetric expression ||] 

This means that 



J^'ip{x = ma) — > —7 



+ (5.14) 



where 7 = fi,^/2ma^, Vm ~ V{ina), and i/jm ~ '4>{raa). In matrix notation this means 
that 



^V = e^, (5.15) 
where the Hamiltonian has the matrix elements 

J^mm' = -lSm,m'±l + (Kn + '^l)5mm' ■ (5.16) 

Interestingly, this is the form of a nearest neighbor tight binding Hamiltonian with 
hopping element 7 and a constant shift by 27 of the on-site energies from Vm to Vm + 
27. For Vm = the solution is the usual cosine band shifted by 27, 
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Fig. 5.2: Plot of the 
parabolic (full line) and 
cosine (dashed line) 
dispersion relations, 
Eqs. 



(5.17 



, |5.18| ). The 
dash-dotted line shows the 
ratio of the two disper- 
sions. For a quarter-filled 
band, e = 7, the ratio is 
9/7r2 ~ 0.91. 



e(fc) = 27(1 - cos fca), (5.17) 
so that e{k = 0) = 0. Expanding the dispersion relation near fc = gives 

£(fc)~7(fca)' = ^, (5.18) 
2m 

and for energies e < 7 the finite differences approximation provides a very accurate 
description of the continuous problem, see Fig. |5.2[ Making the grid finer — a smaller 
— improves the accuracy and/or allows for a treatment of higher energies. 



5.2.5 Retarded Green function 



Since we are considering an open system the Hamiltonian, Eq. (5.16), is in principle 
represented by an infinite matrix. This also means that the retarded Green function 



= [e - ,Jf + iS]-\ (5.19) 

involves the inversion of an infinite matrix. If we just truncate the Hamiltonian ^ 
J^fc with J^fc being an M x M matrix it would correspond to a closed system - the 
wire would no longer be connected to the outside through the leads. To deal with this 
we consider the retarded Green function of the closed system, suppressing the energy 
sub-script. 



(5.20) 
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and the Dyson equation for the full retarded Green function 



(5.21) 



The latter can formally be solved to give 



(5.22) 



which leaves us with the job of obtaining the retarded self-energy E*" due to coupling to 
the leads. Normally this involves approximations but for this particular problem where 
the electrons of each subsystem are non-interacting this is not necessary. To obtain the 
retarded self-energy we write the full retarded Green function as a 3 x 3 block-matrix 



'32 



^13 



'33y 



(5.23) 



where each of the three blocks along the diagonal corresponds to the left lead, the wire, 
and the right lead, respectively. For the Hamiltonian we can of course do the same 



h 

/it JTr, 





,-7 



(5.24) 



where J^i^ and are the Hamiltonians of the semi-infi nite leads and th e coupling of 
the leads to the wire is described by h. Combining Eqs. ( 5.19| , 5.23 , 5.24 ) we get 



^12 
^2'2 



^32 



^13 
^2'3 



^33y 



nr 






(5.25) 



From this result it is straight forward to obtain the sub-matrix ^^22 which is the full 
retarded Green function inside the wire 



'22 



(5.26) 



Comparing to Eq. (5.22) we read of the self-energies S'' = S£ + EJ^ due to coupling 
to the leads 



(5.27) 
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which are given in terms of the diagonal part of the retarded Green functions of the 
semi-infinite leads taken at their respective end-points. 

For a semi-infinite lead (for convenience it starts at site 1) the states are those of a 
semi-infinite tight-binding chain. 



ipi = j= = iV2smkma. (5.28) 

a/2 



By a suitable contour integration it follows that 



and thereby 



{S£}™m' -5i,„JiTO'7CxpiA:a, {I]r}„„„' = -(5A/„JA/™/7CxpiA:a. (5.30) 
Here, k{e) is given by the dispersion relation £{k) in Eq. ( ^.17 ). 



5.2.6 Low-temperature limit 

We will end by summarizing the results for the low-temperature limit where 

A,;H~i-fi^7A„ A, = Re{{<^}^® {<^A<^}}. (5.31) 



Performing the uj integration in Eq. (5.3) we get 



2 



G21 = jC^T^Y {'^siAiC/iaAa}. (5.32) 

This is the main result and it forms the basis for all numerical work reported in 
Refs. [E, F, G, H]. The generalization to several transverse channels is straight forward 
and it provides a very efficient method to investigate drag in different geometries and/or 
for different disorder configurations. 



We note that starting from the Fisher-Lee relation [106| the Landauer conductance 



Gii = dlj/dVi of the individual wires can be expressed in a similar trace-form, see 



e.g. Refs. [B, 113|, 



G,,-^Tr[rL^C.rj,{$#,V}^]. (5-33) 
where the leads are described by 

r; -*[s];-{s;}t]. (5.34) 
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5.3 Two-dimensional formulation 



In the presence of transverse degrees of freedom we can return to Eq. ( 5.7 ) and add the 



transverse integral that was absent in the ID calculation. The net effect is to modify 



the matrix A by a sum over the transverse direction and with the numbering in Fig. 5.1 
we get 



^ Ny (M^-l) + {n-l)M^ 

A^-,,^^,/ ( ]\I iT ^ ^ ^m^rh+l^rhm' <^7nm^m+l,m' • (5.35) 

' ' n=l m=l + (n-l)M^ 



However, this just means that A^^ is an Ny x Ny block-matrix with A in each block 
along the diagonal 



A^^ element: 

, 2D 



A^^ =diag(A,A,A,A...). (5.36) 

The real difference between the one- and two-dimensional systems is thus at the level 
of obtaining the retarded Green function where one need to modify the self-energies, 
see e.g. [0, |t|. Choosing a way of numbering different from that in Fig. ( ^.l| ) does not 
change the overall structure but only that of U12 and A. 



5.4 Implementation 

The matrix formulation is readily implemented on a computer and the accuracy can 
be increased simply by having more lattice points N for a given size of the sample. 
The retarded Green function can be obtained either by the suggestive direct matrix 



inversion, see e.g. Ref. [114|, or by a recursive method, see e.g. Ref. 1115|. The 
numerics to be presented here are based on the former 

For disordered systems the discretized set of equations are often considered a model 
system where the lattice constant is just set to unity and the important pa rameter is the 



band-filling, see e.g. Ref. [ |1 16| |. Alternatively, one can also view Eqs. ( p.32| , p.33[ ) as 
formulas for a tight-binding system where 7 is a 'hopping matrix element' between 
different orbitals and the local potential Vm + 27 is the energy of the orbital localized 
at site TO. 

For the implementation it has been useful to carry out tests for both the Landauer 
conductance Gu and the drag conductance G2i. 

For the Landauer conductance it is obvious to test the numerics for various geome- 
tries which can also be solved analytically, e.g. ballistic regime, single barrier, and 
double barrier tunneling. For disordered systems the outcome of the Anderson model 



can be directly compared to the analytic scaling results of Abrikosov 1 1 17 1 by changing 
the ratio of the mean free path t to the length of the wire L from the localized regime 
{i/L < 1) to the de-localized regime (£/L > 1). 

For the transconductance G21 it is difficult to obtain analytical results. However, 
drag between one-dimensional ballistic wires is an exception where numerics can be 
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compared to semi-analytical results where spatial integrals are evaluated numerically 
by some other method. Another test is the predicted peaking of G21 at the on-set of 
modes in either of the two wires. Especially this can be seen for systems with resonance 
transmission, i.e. peaks in Gu. 



Chapter 6 

One-dimensional wires 



In the following we consider one-dimensional wires which is relavant to the situation 
where the width is comparable to the Fermi wave length so that only a single mode is 
propagating at the Fermi level. We first consider the ballistic regime and afterwards we 
study the fluctuations of Coulomb drag in the quasi-ballistic regime. First we however 
make a unitary transformation by which the current matrix / becomes diagonal. 



6.1 Unitary transformation 

We start from Eq. ( [4.56| ) where in one dimension all matrices / and p are 2 x 2 matri- 
ces. In the presence of time-reversal symmetry {S = S'^) we write the 2x2 unitary 
scattering matrix as 



where TZ = 1— T = |rp = |r'p has been introduced. The current matrix, Eq. ( 4.6C| ) 
can now be written as 



2m 



(t3 - S^T^'S) 



h 

2m 



T 



-r 



Next we employ a unitary transformation 



U = 



(6.2) 



(6.3) 



that satisfies 



m 
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(6.4) 
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by choosing 

\u\' = l{l + VT), 



1 



2V , , 2 

We are free to choose the phases as 

i(0-e)/2 



V = \v\e 



This means that 



(1 - VT), 



u = \u\e 



vu 



-t{4>-e}/2 



(6.5) 



(6.6) 



Tr{/p(x')p(a;")} = -VrTr{r3p(x')p(x")}, p{x) = U p{x)U^ , (6.7) 
m 

where we have inserted 1 = U^U three times and made use of the cycHc invariance of 
the trace. Since the diagonal elements of p are real this also means that 



2h 



ImTr{/p(x')p(x")} = —Vflmp+^ix')p_+{x"). 



From Eq. (4-. 56) it now follows that in the new basis 



A{uj,x',x") = i2Trf h(-^Y huj — VTlmp+^{x')p^+{x"), 
and with the choice of relative phase in Eq. (|5.6|) we have that 



(6.8) 



(6.9) 



n[p^-{x) - p++{x)\ 

+ (1 + VT)p+.ix) - (1 - Vf)p*_,_{x) 



(6.10) 



has no importance in the context of Eq. (5.9) 



though the phase-factor expi((/) - 
In the following we will apply Eq. ( |5.9| ) in our studies of both ballistic and disordered 
wires. 



6.2 Ballistic regime 

In the ballistic regime we have eigenstates given by 

^p^{x) = k-'^/^exp[±ikx], £ = ffk^/2m- p, 



(6.11) 



where the normalization is that of scattering states. In or der to calculate the triangle 
function there are several possibili ties: starting from Eq. ( 5.11 ) we may either i) cal- 
culated / and p and then use Eq. ( 4.37 ) or //) calculate the spectral function £/s and 
use Eq. ( [4.43| ). It can be checked that both routes of course give the same result. Here, 
we will only consider the low-temperature regime and use the approach of scattering 
states, Eq. ( [4-.56 ), as also suggested by the normalization in Eq. (5.11). 
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For the ballistic regime T = 1 and it is thus easy to utilize the unitary transforma- 
tion. Eq. ( 6.1Cl| ) gives, omitting the phase-factor. 



p+-ix) = p*_+ix) = p+-ix) = k 



l,— l i2kFX 



and thereby 



(6.12) 



lmp+_(a;')p-+(a;") = kp^ sin 2k f{x" - x'). 



Using Eqs. (6.9, 6.13) we get 



(6.13) 



A(uj,x',x") = -5-sin2fcF x" - a;'), 



and from Eqs. (4.36, 4.44) it now follows that in the low-temperature hmit 



(6.14) 



G21 — 



1 (fcT) 



h 12 Sp'Sp 



1)^(2) V'^'-F '^i 



(1)^(2) 

^(1)^(2) 



da;'id.4da;'/da;^' 
sin 2fc<J ' (a;'/ - x'^ ) sin 2fc<?> (x^' - Xs). (6.15) 



Here, we have allowed for different Fermi levels of the two wires — the drag is largest 
when they are equal. 

6.2.1 Matching Fermi levels 

For matching Fermi levels we use 



sm2kp(x'[ — x\) sin2fci?(a;2 — 2:2) 



_|_ ^i2kF{x'l~x^) ^-akF{x'^-x'^) _ ^-i2kF{x'^+X2)^i2kF{x[+x'2)l j^g-j 



to introduce the Fourier-like transform 



The drag conductance can then be written as 



G21 — 



e2 1 [kTf |C/-(2fc;^)|'- |t/+(2fc;^)|' 



h 12 4 



2^2 



(6.17) 



(6.18) 



In the following we consider the long-wire limit. 
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6.2.2 Long-wire limit 

For kpL ^ 1 the contribution from [/j^(2fci?) goes to zero due to the rapidly oscillat- 
ing exponential and we get 



eM (fcr)2 \Uu{2kp)\ 
^^^"XT2^| 2^' ^^'^^ 

in terms of the usual Fourier transform of the interaction potential Ui2{q) = Ui2{q)- 

This relation is well-known, see e.g. Ref. and is easy to understand: In 

long translationally invariant wires the electrons can only undergo back-scattering or 
forward-scattering. The back-scattering from k ^ +kp to fc ~ —kp and vice versa 
corresponds to a transferred momentum of q ^ ±2kp to the other layer whereas 
forward-scattering with q ^ does obviously not contribute to drag since no momen- 



tum is transferred. In Fig. 3.1 each of the two interaction lines will carry momentum 
2kp (opposite to each other) and each line contributes a factor Ui2{2kp) to the drag so 
that G21 oc Ul^.i'^kp). 

6.2.3 Long-range Coulomb coupling 

We consider a long-range Coulomb coupling 

[/l2K,X^)=C/([K-.T^)2+d2]l/2^^ [/(^)^^^^ (6.20) 

AneoCrr 

where d is the wire separation. In this case the Fourier transform of the interaction 
becomes 



fL ^iq{x~y) 

I — ^ 

which for long wires, kpL » 1, can be approximated by 



Ui2{q) = U[d) {kpd) kp jj ^ Axdy ; " , === (6.21) 



e 



iqz 



U^2{q)-U{d){kpd){kpL)\ dz^==. (6.22) 



The integral can be performed analytically, see Eq. (3.754.2) in Ref. [118j, with the 
result 

Ui2{q) ~ U{d) (kpL) {2kpd) Ko{2kpd) (6.23) 

where Kq is a modified Bessel function of the second kind of order zero. This means 
that 



G21 ~ f — ^ (kpL) {kpd) Ko{2kpd)] . (6.24) 
n b \ep £p J 
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Fig. 6.1: Plot of G21 for two ballistic wires, Eq. ( |5.18 ), as a function of the length 
kpL for different values of the separation kpd in the case of lo ng-ran ge Coulomb 
coupling. The dashed curves show the asymptotic result, Eq. ( 5.24 ). Note that 

U{kp^) = U{d) kpd. 



Fig. 3.1 shows a numerical evaluation of Eq. (6.18) along with the asymptotic result, 
Eq. (5.24). The long wire limit was also (implicitly) studied in Ref. [[73t], but with no 



justification of the approximation. As seen it gives the right asymptotic behavior, 
but under-estimates the drag conductance by even orders of magnitude for short wires. 
The oscillations for small values of kpL are due to the sin functions in the triangle 
functions, see Eq. ( 5.14 ), and they are not a result of the numerical method that we 
have used. The reason that their magnitude appears so dramatic is due to the log-scale 
plot. 



6.2.4 Screened Coulomb coupling 

In the following we study screened Coulomb interaction where 



f/i2(x'i,4) = tj{[{x'^ - 4)2+^2]'/'), U{r) = t/(r)exp(-r/r,,), (6.25) 

with Ts being the screening length. For d ^ rg drag is strongly suppressed due to 
screening. In the other limit d <C rg we may neglect d in the exponential screening and 
in the long-wire limit get 
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0.5 



Fig. 6.2: Plot of C/i2(g). 
The top curve is for 
long-range interaction, 
Eq. (5.23), and the two 



other are for screened in- 
teraction, Eq. ( |5.26| ), with 
rs/d = 10 and r^/d = 5 
from above. For rg oo 
Eq. ( |5.26 ) approaches the 
result of Eq. ( |5.23 ). 



Ui2{q) ~U{d){kFL){kFd) I Az^ 



^/z^^ + d^ 

= U{d) [kpL) {kFd)7rRc{Ho{iqd + d/rs) -Yo{iqd + d/rs)}, (6.26) 

where _ff„ is the Struve function of order n and F„ is the Bessel function of the second 
kind of order n. This means that 



^e^ 1 kT U{d) 

XTiRc{Hi:,{i2kFd + d/rs)-YQ{i2kFd + d/rs)]\ . (6.27) 



In Fig. 5.2 we compare Ui2{q) to the case of long-range interaction. For moderate 
screening (r^ > d) the effect of screening is mainly to modify the short wave length 
behavior and since drag probes U12 {^kp) the effect of screening is not dramatic. How- 
ever, as we shall see later the fluctuations of drag probes U12 {q ^ 0) and then screening 
will make a more pronounced difference. 



6.2.5 Short-range Coulomb coupling 

We consider a point-like coupling 
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Ui2{xi,X2) = 'W{d)6lkp{xi - X2)]- 

This gives 

Ui2iq)^kFL^id), 

SO that in the long wire limit kpL ^ 1 

1 1 fkT 



G21 — 



■ kpL 



(6.28) 



(6.29) 



(6.30) 



h 122 yep ep 

All three models for t he interaction ha ve the common feature that G21 oc in the 
long-wire limit, see Eqs. (6.24,5.27, 5.3C). 



6.3 Quasi-ballistic regime — analytics 

We consider wires with a mean free path t much longer than than the length of the 
wires, ^ 2> so that they behave as almost ideal quantum wires with the Landauer 
conductance of the two wires characterized by 



{Gu) ^Gu{l^oo) = {{6Guf)^Q, ^ » L. (6.31) 

Though the mesoscopic fluctuations of the Landauer conductance Gu are vanishing in 
the quasi-ballistic regime L ^ £ it turns out that the fluctuations of the drag conduc- 
tance G21 are finite and they can even exceed the mean value 



(Gai) ^G2i(^^oo), {{5G2iff'^ ^ {G21), ^ » L. (6.32) 

In the following sections we derive this relation and study it both perturbatively and 
numerically. 

We consider long wires, kpL ^ 1, in the low-temperature regime. We are going 
to consider the magnitude of the fluctuations relative to the mean value and thus we for 
simpUcity use the short-hand notation 



G2i=C / U i2{x, y)U i2{x' y' ) l^i{x, x ) l^2{y' ,y), 



(6.33) 



with, see Eq. (6.9), 



l\{x,x') = l-ai{~p+_{x)~p.+ {x')]. 



(6.34) 



Here, / means an integral over all spatial degrees of freedom — in this case x, x' , y, 
and y'. The frequency integration has been performed to give the T^-dependence and 
all prefactors are included in C, i.e. C cx [e? / h){kTY . 
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6.3.1 Ballistic regime 



The notation is most easily illustrated by re-considering the ballistic regime. Since we 

— 1/2 

only have a single sub-band we can safely incorporate the kp normalization, see 
Eq. (6.11), into the prefactor C. With this notation we get 



A(a;,a;') = Imjf 



} = &\n2kF{x' — x), 



(6.35) 



which can be compared to Eq. ( |6.13 ). The over-line of the triangle function is used to 
indicate the ballistic regime. Eq. (6.19) now becomes 



G2i(oo) = lim G21 W = C-UU2kp), 



with 



dxidx2 e'«("^""=^C/i2(a;i,a;2), 



being the Fourier transform of the interaction. 



(6.36) 



(6.37) 



6.3.2 Mutually un-correlated disorder 

In the presence of disorder we consider the correction to the triangle function. 



A{x,x') ^ A{x,x') + 5A{x,x'), (6.38) 

and if the disorder potentials of the two wires are mutually un-correlated (uc) we can 
write 

(Ai(x,a;')A2(2/,y'))ue = (Ai(x, x'))(A2(2/, y')), (6-39) 

as 

(Ai(.T,x')A2(2/,y')),, = Ai(x,x')A2(y,y') + Ai(x,x')(^A2(y,y')) 

+ {SAiix,x'))A2iy,y') + {6Ai{x,x')){6A2iy,y')). (6.40) 

In the quasi-ballistic regime we can make an expansion in the disorder strength. Dia- 
grammatically we have that the average of a triangle function with only a single impu- 
rity line attached vanishes during disorder averaging since for the averages over disor- 
der we assume a potential with zero mean and a short-range correlation 

{V,ix)) = 0, {V,{x)V,{x'))^^ = S,,V^d{x - x'), £ = [hvFlV^f. (6.41) 

Here, the back-scattering mean free path £ is related to the disorder strength Vb within 
the Born approximation. This result is easy to show with the aid of Fermi's golden rule 
and we will do this calculation in the discrete case at a later point. 

This means that all corrections to averages (A.;) and (A; A.;) will be of second or- 
der in the disorder strength, i.e. to lowest order in 1 /kpl. Diagrammatically this means 
that there should be impurity lines connected to two points on the triangle function(s). 



6.3. QUASI-BALLISTIC REGIME — ANALYTICS 67 




Fig. 6.3: Examples of the lowest order connected diagrams in the diagrammatic ex- 

1/2 

pansion for the fluctuations giving rise to a ( ((5G2i )^ ) oc Vq oc l/kp£ dependence. 
Due to momentum conservation diagrams with only one impurity line do not con- 
tribute significantly to the drag fluctuations since they have to carry momentum q ~ Q 
corresponding to forward-scattering. The diagrams (a) and (c) are relevant for both 
correlated and un-correlated disorder whereas the diagrams (b) and (d) are relevant for 
correlated disorder only. 



6.3.3 Mean value 

For the mean value we note that 

{5A{x, x')) < A{x, x'), L < £, (6.42) 
so that the corrections to the mean value are small, i.e. 



(G2i)^^ ~ G2i(c^) + oii/kpe) = c-u^^{2kF) + oii/kpe). (6.43) 

To lowest order in 1/kpt the effect of disorder is just to broaden the distribution of 
the drag conductance without shifting it. We now calculate the magnitude of these 
fluctuations as a function of the mean free path. 



6.3.4 Fluctuations 

For the fluctuations, 



{{5G2if) - {Gl,) - (G21)', (6.44) 

the presence of even weak disorder makes a big difference because the broken transla- 
tion invariance allows for transferred momentum different from 2kF. 

Fig. shows different examples of diagrams with two impurity lines. The im- 
purity lines carry momentum 2kF corresponding to disorder induced back-scattering 
within the wires. Diagrams with a single impurity line are in principle also possible. 



but only if the line carries momentum q = Q corresponding to forward-scattering [ 1 1 9 1 



However, forward-scattering in the disorder channel will not contribute significantly to 
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the fluctuations for a sufficiently long-range interaction and for a short-range inter- 
action where forward-scattering becomes important the fluctuations are substantially 
suppressed when compared to the mean value. 

The main difference between diagram (a) and (c) is that in (a) the interaction lines in 
each sub-diagram must carry the same momentum whereas in (c) one of them can carry 
e.g. zero momentum while the other carries 2fci?. Each impurity line is characterized 
by ~ L/^ and the interaction lines by C/i2((7). For the fluctuations the contribution 
from diagrams of type (a) can be estimated as 

(a): (((5G2l)2)^,^^C2(7^l)Vl42(2AF), (6.45) 
whereas for diagrams of type (c) we get 

(c): {{SG2if)^^--C^ni){n2)Uf2{2kF)U?2{0)- (6.46) 

Since Ui2{0) ^ C/i2(2fci?) for a long-range interaction this means that diagrams of 
type (c) give the major contribution to the fluctuations. This can be tested numerically 
by noting that the diagrams of type (a) are relevant to the case where wire 1 is disor- 
dered and wire 2 is either disordered or ballistic whereas diagrams of type (c) are only 
relevant to the case of both wires being disordered. Indeed, by numerically calculating 
the fluctuations for a system where one of the wires is ballistic and the other disordered 
we have found a very dramatic reduction of the fluctuations compared to the case of 
both wires being disordered. 

To show the result in Eq. ( p6| ) in more detail we note that in the quasi-ballistic 
regime 



{{6G2if)^^ c^C^ J Ui2{x, y)Ui2{x', y')Ui2ix, yW^x' , y') 

X (5Ai(x,x')^Ai(S,a;'))('5A2(y',y)^A2(r,y)), (6.47) 

corresponding to diagram (c) in Fig. |6.3| where all four triangle functions have an im- 
purity line attached, i.e. two points for each of the averages. 

6.3.5 Lippmann-Schwinger equation 

The task is thus to calculate the correlation ((5Ai(.T, x')(5Ai(a;, x')) to lowest order 
(second order) in the disorder strength. For quasi-ballistic wires, where T = 1—7?. ~ 1, 
Eqs. (plpO|)give 



A(x,x') « Im{p+_(x)p_+(x')}, (6.48) 
and formally we thus have that 

(5A(x, x') = Im{5p+_(a;) p^+{x') + p+-{x) Sp^+{x')] 

= Im{d>+_(.T)e*2fe^^' + e-'^'"=-'=SpX_{x')}. (6.49) 
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Using the Lippmann-Schwinger equation, see e.g. Ref. [ |12C| , we have to lowest 
order in the disorder strength that 



where 



Jo 



(6.50) 



(6.51) 



is the unperturbed retarded Green function, i.e. that of a ballistic wire. Introducing 

(x) = u(x) + v{x) this means that 



u(x) 
v(x) 



1 



ihvp 
1 



ihvp 



dx^x)- / AxVix) 



The disorder potential V{x) is "self-averaging" in the sence that 

dxV{x) - {a-b){V) = 
and thus the last term gives a small contribution 

6A^-'\x,x') ^ 



(6.52a) 
(6.52b) 

(6.53) 



ihvp 



+ dxVix)- J ^ dxVix)jA{x,x'), 



which we will neglect so that 5p^ (x) u{x). 

For long wires, kpL ^ 1, we neglect terms oscillating with ikp and get 



{u{x)u{x)) ^0, (u{x)u*ix)) ^ (^1 



2\x~x\ 
L 



(6.54) 



The prefactor L{Vo/hvF)'^ = L/l can be related to the reflection coefficient of the 
wire. To see this we consider the Dyson equation for the retarded Green function to 
second order in the disorder 



{x,x')^%{x,x')+ / dx%{x,x)V{x)%{x.x) 



dxdx'%{x,x)V{x)%{x,^)V{x!)%U.x'), (6.55) 
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and from the Fisher-Lee relation 1 106| we notice that 



n = i-T = 1- {rwFfp''{L, 0)1 



It then follows that to second order in the disorder 

2 



\hvF ) i 



(6.56) 



(6.57) 



This means that 



(5 A(2;, x')i5A(x, x')) w (7e){/(;r' - i', a; - x) + /(x - x, x' ~ x') 

- f{x' -x,x-x')- fix - x', x' - x)}/2, (6.58) 



where we have introduced the function 



/(z, z') = (1 - 2\z\/L) cos2kF{z'). 



(6.59) 



For the fluctuations we thus get 4^ — 16 terms, but only 4 of these terms wit h the form 
cos 2kF{x' — x') cos 2kF{y' — y') are finite for long wires, fci?i ^ 1. Eq. ( 6.47 ) now 
becomes 



((JG2i)2) ~c2(7^l)(7^2) / c/i2(x,2/)t/i2(x',j/')t/i2(s,y)t/i2(.T',r) 



X f{x - X, x' - x')f{y -y,y' - y'), (6.60) 



and introducing 



we get 



d.Tida;2dx2da;2 



((-^621)2)^^ c'-{n,){n,)uU2kF)uUo)- 



(6.61) 



(6.62) 



This is in close agreement with the form suggested in Eq. ( |5.46 ). To estimate the 
relative magnitude we divide Eq. ( 5.62| ) by Eq. (5.43 ) and get 



(('5G2i)2);f [2{n,){n2)uh{2kF)u?, 



I 1/2 



(G21), 



Ut^{2ki 



(6.63) 
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Fig. 6.4: Comparison of J7i2(0) to J7i2(0) for a long-range interaction. The plot shows 
a numerical evaluation (indicated by the symbols) of the ratio ?7i2(0)/?7i2(0) as a 
function of kpL for three different separations kpd — 1, 2, 3. The dashed lines are 
guides to the eye. 



From Eqs. ( 5.43| , 6.63 ) it follows that the relative magnitude of the fluctuations is of the 
order 

(G2l)^^ Ul2[2kF) 



since [/i2(0) ^ f7i2(0) within a numerical factor of order 2, see Fig. 5.4. This is the 
key result on one-dimensional wires reported the first time in Ref. [E] and later studied 
in more detail in Refs. [F, G, H]. 

Even though (TZ') <C 1 the effect of long-range interaction may compensate this 
since J7i2(0)/C/i2(2fcir) 3> 1, see Fig. |6.2| . This results in relative fluctuations of order 
unity, i.e. fluctuations comparable to the mean value. As we shall see in the next 
section, Eq. ( |6.63| ) fully accounts for our numerical results. 

6.3.6 Mutually correlated versus mutually un-correlated disorder 



We will now relax the assumption of mutually un-correlated disorder, see Eq. (6.40). 
In fact we will study the other extreme limit where the disorder potentials of the two 
wires are mutually fully correlated (c), i.e. 



Vi{x) = V2{x), {v,{x))=Q, {V,{x)Vj{x'))^^Vo^6{x~x'). (6.65) 
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Fig. 6.5: Diagrams for the mean drag with correlated disorder, (G21) ^- The first term 
is the ballistic limit and the second term shows an example of the lowest order diagram 
in the diagrammatic expansion of the mean value with a.V^ oal jkpl dependence. 



The calculation of the corrections to the mean value in the case of fully correlated 
disorder is not that different from the calculation of the fluctuations in the case of 



un-correlated disorder. Fig. p.5| shows an example of a lowest order diagram which 
contributes in the case of correlated disorder Again, we consider back scattering in the 
disorder channel so that in the second term the impurity line carries momentum 2kp. 
This allows one of the interaction lines to carrying also 2kp while the other one carries 
zero momentum, g ~ 0. For the mean value we thus get the estimate 



(G2l)^ ~ Ul2{2kp) + (7^) C/i2(0)t/i2(2fc;^), 

and from Eq. ( |5.43 ) it thus follows that 

(G2i)^ ~ [1 + (7^) (7i2(0)/[/i2(2fcj.)l X (G21) 



(6.66) 



(6.67) 



Since (7?.) C/i2(0)/?7i2(2fcF) can be of order unity for long-range coupling it means 
that correlation of the disorder may enhance the mean value by up to a factor of order 
two compared to un-correlated disorder. Similarly to the case of fluctuations this effect 
is due to a combination of back scattering in the disorder channel and forward scattering 
in the Coulomb channel. The case of correlated disorder was recently studied for bi- 



layer systems where it was also found to enhance the mean value [ 104| 



For the fluctuations the calculation can be carried out using simple diagrammatic 



arguments. We again consider the diagrams in Fig. 6.3. For correlated disorder both 



of the diagrams (c) and (d) contribute equally whereas for un-correlated disorder only 
the diagram (c) is relevant. More generally, for each topologically different diagram 
contributing in the case of un-correlated disorder there are two similar diagrams con- 
tributing equally in case of correlated disorder Of course there are also other possible 
diagrams in case of correlated disorder but to lowest order in l/kp£ this means that 



--{a) 



{{SG2^r)^ ^{a)+{b) +{c)+{d) 



((5G2l)2)^ 



(a) + (c) 



(6.68) 



where (a), (6), (c), and (d) refer symboUcally to the diagrams in Fig. |6.3| averaged 
over disorder This means that the fluctuations in the case of correlated disorder will 
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be enhanced by a factor of >/2 compared to the case of un-correlated disorder, 

((JG2i)2)f ^V2x((5G2if);f. (6.69) 

Though these arguments may at this point seem simpUstic, the result turns out to be 
correct as we shall see later when comparing to numerics. 



6.4 Quasi-ballistic regime — numerics 



The study of the statistical properties of disordered systems is an interesting example 
where Eq. ( 5.32 ) can be applied to a large ensemble of different disorder configurations. 
For the di sorder we n eed a way to represent a short-range correlated disorder potential, 
see Eqs. ( |64l||6!65| ). Here, we use the most often used model; the Anderson model 
[121] which the first time was introduced in the study of one-dimensional localization. 



6.4.1 Anderson model 



In the Anderson model with diagonal disorder [121, 116] the site potentials Vni at site 
n of a given wire i are statistically independent with each site energy taken from a 
uniform distribution of width W centered around zero, i.e. 



p{Vm) = e{W/2 ~ Vn,)/W. (6.70) 

This means that 

(14,:) =0, {VmVn'j)^=S,,n'-^, {VmVn'j)^^^5:j{VmVn'j)^. (6.71) 

Again, the disorder strength W can be related to the back-scattering mean free path. 
For a tight-binding chain with N sites and no disorder we have eigenstates V'^j = 
no). We calculate the rate for back-scattering from Fermi's golden rule, 
suppressing the wire index i. 



(6.72) 



where the disorder potential V is treated as a perturbation. The factor following the 
delta function restricts the k' sum to those relevant for back-scattering, i.e. k' opposite 
to k. For a chain with large N we take the sum into an integral and since 5{e — e') — 

\hvk\~^ [S{k - k') + 5{k + k')] we get 



r(fc) Nh^\vk\ 



(6.73) 



The corresponding mean free path is given by t{k) ^ Vk/{T ^ (fc)) and for the disper- 
sion in Eq. (5.17) we get 
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£=12a{'ijSF-sl)/W'^, (6.74) 

at the Fermi level. This result agrees with Ref. jl 16[ ] except for the trivial constant shift 
of the energy by 2j. 

6.4.2 Numerical results 

We will first consider two special cases: i) both wires being disordered, but with Vi 
and V2 mutually un-correlated (as in Ref. [E]) and ii) both wires being disordered and 



fully correlated, i.e. Vi = V2 (as first suggested in Ref. [ 104 |). Finally, Hi) we consider 
the more general case where the two wires have partly a correlated (common) disorder 
potential 14 and partly an un-correlated Vuc one. 



Vi=V, + Vuc.i, l^2 = Ve + Kc,2, W^ = W^+Wl, (6.75) 

but with the impurity strength adjusted so that the mean free path is the same as in the 
two special cases (the last equation). 

We consider quarter-filled bands {ep =7) and wires with = 100 lattice points 
so that kpL = (tt/S) x 100. The separation is kpd = 1 and we assume long-range 
Coulomb interaction. 



47reoer\/(? 



(6.76) 



For the case £ ~ 36L (this corresponds to W ^ ep/10) the disorder has as ex- 
pected [117| almost no effect on the Landauer conductance, i.e. (Gu) ~ 2e^//i with 



vanishing fluctuations. However, for the transconductance the situation is very differ- 
ent. Panel (a) of Fig. 6.6 shows a typical histogram of G21 {£) / G21 (00) for I ~ 36L in 
the case of un-correlated disorder. Depending on the disorder configuration G21 {£) can 
be either higher or lower than in the ballistic regime. We emphasize that the histogram 
peaks close to the ballistic value and not at zero drag, i.e. the mean drag is finite and 
positive in agreement with Eq. ( |5.43 ). The variance is of the same order as the mean 
value so that sign reversal for some disorder realizations is possible. The later is repre- 
sented by the negative tail in the histogram. The sign of the drag is thus arbitrary in the 
sense that both positive and negative drag can be observed. 

For the distribution of G21 in the case of un-correlated disorder we find that 



a ~ 1.4, (6.77) 

fits surprisingly well to the data. In fact we have performed these fits to histograms for 
kp£ ranging from 10"^ to 10^ with kpL in the range from 100 to 300. By rescaling of 

G21, 



Puc(a;) (X exp 





X — Xq 




Xl 
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4 5 6 7 




Fig. 6.6; Histograms for 
G21 {£) normalized by the 
ballistic result 6*21(00). 
The histograms are based 
on > 10^ disorder config- 
urations and for all three 
histograms the mean free 
path is ^ = 36i. Panel (a) 
is for the situation of mu- 
tually un-correlated disor- 
der, panel (b) for partly 
correlated disorder (Wc = 
Wuc), and panel (c) is for 
mutually fully correlated 
disorder. 



1 01234567 

G2l(^)/G'2l(00) 




Fig. 6.7: Rescaling of the distribution of the drag conductance G21 for mutually un- 
corrected disorder Distributions for kpL ~ 100 and kpd ~ 1 for different values of 
the mean free path kpl- 



G 



21 



G21 — {G21) 
{{5G2iY) 



1/2-' 



(6.78) 



it is possible to let all histograms fall onto the same curve. There are three parameters 
which may be varied: L, and d. For simplicit y we illustrate this by fixing two of the 
parameters and varying the third one. Figs. 5.7, 6.8, 5.9 show the rescahng for fixed L 
and d, fixed i and d, and fixed £ and L, respectively. 

Note that the distribution is non-universal in the sense that it depends on the range 
of the interaction Ui2- It should also be stressed that the particular form i n Eg . ( 5.77 ) 
is not unique and that other similar functions might be fitted equally well 



119|. 



In Fig. 5.1c w e com pare the results of the Anderson model for the case of a uniform 
distribution, Eq. ( |5.7C| ), to a Gaussian distribution, but with the same mean value and 
variance. As seen the distribution of the drag conductance G21 does not depend on the 
particular choice of disorder model. 




Fig. 6.8: Rescaling of the distribution of the drag conductance G21 for mutually un- 
corrected disorder Distributions for different values of kpL for kpd = 1 and a fixed 
ratio of kpi/kpL. 



For the same system parameters but now with fully correlated disorder we get a 
very different distribution as seen in panel (c) of Fig. 6.6 Again we find that rescaling 
is possible, see Fig. 6.1 1 



In this situation the mean transconductance is enhanced by approximately a factor 
of two compared to th e un-correlated case. This confirms, at least qualitatively, the 
prediction of Eq. ( 6.67 ). The fluctuations are also enhanced. However, since the mean 
transconductance is enhanced by almost a factor of two compared to the ballistic limit 
this also means that there is no disorder configurations giving rise to negative drag. 

In Fig. 6.12 we show the dependence of the fluctuations on the mean free path 
£ which has the expected l/kpt dependence. We also notice that correlated disorder 
gives rise to slightly larger fluctuations compared to un-correl ated d isorder In fact they 
are exactly enhanced by a factor of V2 as predicted, see Eq. ( 5.69| ). 

For partly correlated disorder we as expected get fluctuations which are enhanced 
by a numerical factor between unity and \/2. For the histograms we observe how 
the, say, histogram in the case of un-correlated disorder evolves into that of correlated 




Fig. 6.9: Rescaling of the distribution of the drag conductance G21 for mutually un- 
corrected disorder. Distributions for different values of kpcl for kpi = 36 x kpL and 
kpL - 100. 



disorder when increasing the degree of correlation, see panel (b) of Fig. 6.6. It should 
however be stressed that for partly correlated disorder the histogram is not a simple 
superposition of the histograms in the two extreme limits. 

To summarize we find that our perturbative results for the fluctuations are in ex- 
cellent agreement with our numerical studies. Furthermore, the numerical studies give 
insight into, not only the statistical moments of G21, but in fact the full distribution 
P(G2i) which we presently have no analytical predictions for 
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Fig. 6.10: Rescaling of the distribution of the drag conductance G21 for uniformly and 
Gaussianly distributed on-site energies. 
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Fig. 6.1 1: Rescaling of the distribution of the drag conductance G21 for mutually cor- 
related disorder. Distributions for kpd — 1 and kpL ^ 100 for different values of 
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Fig. 6.12: Plot of fluctuations {[5G2i{e)]'^) as a function of the mean free path kpi. 
In both cases the expected 1/kpi behavior is born out by the numerical calculations. 
The full lines are Eq. (5.63) with no free parameters. The expected enhancement. 



Eq. ( 5.69 ), for correlated disorder by a factor of V2 compared to un-correlated disorder 
is also confirmed by the numerical calculations. 



Chapter 7 

Quantum dots 



A quantum dot is a sub-micron cavity with electronic properties governed by the phase- 
coherence of quantum mechanics. It can e.g. be one of Nature's own small metallic 
nano-particles or it can be man-made by confining the electronic degrees of freedom in 
a semiconductor Quantum dots exhibit many exciting properties both in their charge 
transport and in their optical response. Quantum dots with the classical motion being 
chaotic is also a field of tremendous interest often referred to as quantum chaos. 

Here, we consider two Coulomb coupled semiconductor quantum dots such as sug- 
gested in Fig. L5. Our aim is to calculate the statistical properties of the drag conduc- 
tance with the aid of random matrix theory. 



7.1 Statistical averages 

Reviews of the statistical theory of disordered and chaotic systems are given in Refs. 



122, 123 1 . For the appUcation of random matrix theory to the current matrix we will 



mainly follow the notation used in the review of Beenakker | 



7.1.1 Mutually un-correlated disorder 



We consider an ensemble of mesoscopic chaotic systems, as suggested in Fig. 1.5 
We assume that the region where the subsystems couple by Coulomb interaction have 
mutually un-correlated disorder so that, suppressing the integration variables. 



(G21) (X J C/i2C/i2(Ai)(A2), (7.1a) 
(G^i) cx J C/i2C/i2C/i2C/i2(AiAi)(A2A2). (7.1b) 



We consider the low temperature limit and start from Eq. (4.46) — at some later point 
we turn to the formulation in Eq. ( 4-.56| ). The task is to calculate (A;) and ^A^Ai) 



which we do using random matrix theory. 
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7.1.2 Correlation of current and wave functions? 



To lowest order in l/kpi the average of two wave functions is, see e.g. Ref. [124], 

{(t)*^{x)(t)s{y)) = 6js{(f>yix)(l>yiy)) « 5^,5^ (|27r^ (/)*(a;)(?!)^(y)5((£y)^, (7.2) 

7 

where Se is the mean level spacing. Introducing the spectral function 

^(r,r') = 27r^C(r)</)„(r')<5K), (7.3) 

a 

we get 

{^;{x)^siy)) ~ S^s^{^{x ~ y)), (7.4) 

where the average spectral function {£/{x,y)) = — y)) at the Fermi level is 

given by, suppressing the energy subscript, 

{£/{r)) ~ (TO/2ft2)cxp(-r/2€)Jo(fci.r). (7.5) 

Next, consider the average 

{Iai3(l)*{x)(f>s{y)) cx {d:^^ - d:^^){(l)l{xi)(l)p{x2)(f>*{x)(f>s{y)) 

where to lowest order in l/kp£ 

{(I)*a{xi)(l)l3{x2)(p*{x)(j)s{y)) ^ {(f)l{xi)(f>i3{x2)){(p*{x)(l)s{y)) 

+ {raix^)q^s{y)){M^2)r-,ix)). (7.7) 

Due to current conservation the points xi = X2 can be chosen freely. Taking them 
to be outside the chaotic region, the decay of the spectral function makes the second 



(7.6) 

Xi=X2 



term in Eq. ( [7.7| ) vanish if x or y are inside. The first term amounts to performing the 
average over lap and 4i*{x)(j)s{y) separately. Similar arguments hold for higher order 
averages so that to lowest order in 1/kpl 



^ ' Q/37 



(7.8) 



and since / and p are Hermitian 

{A,{iu,r,r')A,{Q,s,s')) 



huj tub {2TTYh? (2i) 
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For the averages we have used that eigenvalues and eigenfunctions are un-correlated so 
that there is quite some freedom in where to place the delta functions. We note that in 
these expressions the triangle functions are only given to lowest order in uj so that the 
right hand sides do not depend on uj. 



7.2 Random matrix theory 



In order to calculate (la-y) and {layla^y) that enter Eqs. (7^, 7^) we apply standard 
random matrix theory for chaotic quantum dots By chaotic we mean that the 
classical motion in the dot can be regarded as chaotic on larger time scales than the 
ergodic time; the time needed for a classical trajectory to explore the entire phase- 
space. We assume that the dot is coupled to the leads via ballistic point contacts with 



N modes. Then the distribution of the 2N x 2N scattering matrix S, see Eq. ( 4. 50 ), is 
given by the circular ensemble; the scattering matrix S is uniformly distributed in the 
unitary group and is only subject to the constraint of a possible time-reversal symmetry. 



7.2.1 Circular ensemble 

Now we consider a unitary M x M matrix U (UU^ = 1) which is distributed uniformly 
over the unitary ensemble with no further constraints; for M = 1 this is to say that 
U = exp i9 is uniformly distributed over the unit circle. Here, we will only need a 
few results which can be found in Refs. |125| ]. The most important one is that if we 
consider a polynomial function 



fiu) = f[Uo.^a^f[u;^,^, (7.10) 

i—n m—1 

then the average (/(/7)) is zero unless p = q and the sets of "left" indices coincide, 
{«„} = {/?n}, and similarly for the "right" indices, {a„} = The first constraint 

means that 



(C/) = 0, (C/2) = 0. (7.11) 
The other constraints are most easily illustrated for p = q = 1 where 



For p = q = 2 the result is slightly more complicated 



(7.12) 
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(7.13) 



General results for p = g = 3 and p = (7 = 4 have been given by Mello [125|. 
We will need the expression for p ~ q ^ A but since it contains (4!)^ = 576 terms 
its application is not a simple task. Instead we consider the Gaussian approximation 
which is the leading-order term in the average 



{I{U)) ^ M-^^.Y.W^'-.PvuA.bvu) +0{M-P-'), (7.14) 

V j=l 

where V is all the pi permutations of the numbers 1, 2, ... p. 

7.2.2 Current matrix elements 

The task is to calculate (la-y) and (la-yla^) but since we, as in Eq. I^T^, may sum 
over the indices and divide by the number of modes we might as well consider 



aj 0:70:7 

Starting from the current matrix 

'-i-Jr'-'^VS]. ,7.16, 
the first average becomes 

07 07 

and similarly it can be shown that 

(^-7^57) = (^) E {{S^r^S}o.^{S^T^S}ay). (7.18) 

a7Q7 Q7Q7 

We will first consider the case of time-reversal symmetry (usually denoted by the 
symmetry index f3 = 1 [|]|) and then afterwards discuss the effect of breaking time- 
reversal symmetry by a weak magnetic field (/3 = 2). Finally, we summarize the results 
in the asymptotic limit ^ 1. 
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7.2.3 Time-reversal symmetry 

In the presence of time reversal symmetry S = which motivates the introduction 



of a unitary 2N x 2N matrix U by S = UU . From Eq. ( 7. 13) we get 



vxy 



X V 



where AI = 2N. This means that Eq. (7.17) becomes 



i^'^l) B=l ^ 



a7 1 



(7.20) 



Similarly, 

= E {U^yUiyU,uU^uU:,Ul,U;M:) (7.21) 

xyzu 



The average is zero unless for the indices 

{yyuu\ = {xxzz}, {vjv^} ~ {vavd}, 



(7.22) 



which in principle allows for simplifications. Here we will just consider the Gaus- 
sian approximation, Eq. (7.14). The calculation is in principle simple although long. 
However, with the aid of Mathematica Eq. (7.18) becomes 



h V M{M'^ + 3AP + 6M) 



2m 



AP 



(7.23) 



The same calculation based on the full expression of Mello [ 125 1 is of course possible 



as soon as the (4!)^ = 576 terms have been fed into the computer. The 4! = 24 terms 
in the Gaussian approximation are on the other hand easily generated by the computer 



itself, see Eq. (7.14) 



7.2.4 Breaking of time-reversal symmetry 

In the absence of time reversal symmetry S is unitary with no further constraints so 
that Eq. (7.12) gives 
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(S^T^S)^^ = J2 rLiU.pU:^) = ^S^p Trr^ = 0, (7.24) 



and thus Eq. ( P.17[ ) gives 



(7.25) 



In the same way 

[5^a5^a + S^oS-ya^^a] , (7.26) 



M(M2 - 1) 



which means that 



\ L J a7 L J a^y / 

From Eq. ( [7.18[ ) it now follows that 



1 



il/2 - 1 



1 



M(A/2 - 1) 



/3=2 V2m/ A'/ + l 



' Af2 



(7.27) 



(7.28) 



Using the Gaussian approximation, Eq. (7.14), instead we get the result (?i/2m) M 



7.2.5 Asymptotic limit 

From Eqs. ( 7.20| , 7.25 ) we have 



p=i,2 ~ 2m'^"''' ^ (-^"7)^=1,2 ~ ^ 



(7.29) 



Q7 



independently of Af = 2N and /3. In fact, the last equation holds not only in the 
asymptotic limit M ':$> 1 but also for any arbitrary number of modes. From Eqs. 
(^^BlOsbwe get 



which also does not depend on /?. From Eq. ( [7.29[ ) it now also follows that 

(/a7^o7>/3=i,2 = const. X r^^T?^ + {{SW^S}ay{S^r^S}a^)^^^^^ (7.31) 
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which in the asymptotic Hmit reduces to 

A'/ > 1. (7.32) 



This is our main resuh of random matrix theory and the result does not depend on (3. 
For the fluctuations that we are going to calculate this means that there will be no effect 
of breaking time-reversal symmetry by a weak magnetic field in contrast to the case of 
universal conductance fluctuations JTl] , [T^ where there is a factor-of-two reduction of 
the fluctuations when breaking time-reversal symmetry [^]. We notice the misprint in 



Ref. [E] where the factor M ^ in Eq. ( 7.32 ) has a wrong power which unfortunately 
propagates throughout the paper. 



7.3 Statistical properties of the triangle function 



7.3.1 Zero mean value 



Sinc e (i g-y) oc t^^ = r^^^Sa^ we can in the second average in (A) replace 7 by a, see 
Eq. (7.8). The second average does not depend on the direction of the current and thus 
we get (A) = which in turn gives 



(G2 



0. 



(7.33) 



This is the first out of two key results for quantum dots. Experimentally, ensemble av- 
eraging could be performed similarly to the case of universal conductance fluctuations 
by magneto-conductance exp erim ents [126|, by varying the shape of the dots [127], or 
by changing the Fermi level 1 128]. The prediction of zero average drag can be consid- 
ered as a test of the degree of ergodicity. Our second key result concerns the finiteness 
of the fluctuations which we study below. 



7.3.2 Finite fluctuations 



Let us now calculate the fluctuations by using Eq. ( 7.32 ) in Eq. (7.9). The first term 
of Eq. ( 7.32 ) does not contribute for the same reasons as for the m ean value. For the 
second term we perform the 7 and a sums according to Eq. ( 4.55 ) and since the sum 
and the energy integral in Eq. ( 4.55 ) can be performed independently due to eigenstates 
and eigenvalue being un-correlated we get 



a/37(i/37 a/37 



(7.34) 



With the Kronecker deltas Sa^Sa-y in the second term of Eq. (7.32) we can now perform 
the 7 and a sums which gives 
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{A,{cj,r,r')A,{Cj,s,s')) 1 

ftcj ftcZ- (277)4^2 - 4(27r)2;i2(2i)2(2iV)(27r)4 

x|(£/,(r,s')M(r',r)M(s,r-')M(s',s)) 

-(M(r-, s)£/,{r', r)M(s, s')M(s', r')) 
-{£/^{r, r')£/^{r', s')^i(s, r)M(s', s)) 

+(^,(r,r'K,(r',sK,(s,s')^,(s',r))}. (7.35) 
Here, we have introduced the spectral function at the Fermi level . 



7.4 Fluctuations of drag conductance 

To lowest order in l/kpi we replace each spectral function by its average and use 
{s^{r, r')) = (j2/(r', r)), so that 



TT'^hLohQ ~ 2(2iV)(27r)4 ' ^ - ' 

where 



J-,(r,r',s,s') = (M(r,r'))(M(s,s')) 



Since 



=0 



(7.37) 



(7.38) 



we perform the oj-integration in Eq. ( 4.36 ) and get 



((5G2if)'/'c.^i^[y't/i2(ri,r2)[/i2K,r^)C/i2(si,^2)C/i2K,4) 



1/2 



(7.39) 



To obtain an estimate of the fluctuations we note that jz/ in Eq. (7.5) is largest for 
kpr < 1 and on long length scales 



Using that approximation for all spectral functions is too crude, but 



(7.40) 
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X [S{r - s)S{r' - s') - 5{r - s')6{r' - s)] , (7.41) 

still gives a finite answer 



1 



h 2,2»N {kpiy \eF 



kT 



[Ur2{ri,S2)Ui2{r\,s'^)-Ui2{r[,S2)Ui2{ri,s'^)y 



1/2 



(7.42) 



Due to the peaked behavior of the spectral functions the mixed terms give a vanish- 
ing contribution and introducing the new variables 



we get 



V ^ S2- S2: 



kT 



(7.43) 



h 32»N (kpe)^ \eF 



2A I {j^{u,0)y{^{v,0)yU^^{u + y,v)U^^{y,0) 



1/2 



, (7.44) 



where A is the interaction area. Assuming that the screening length <C £, the range 
of is longer than the range of [/12 and we further approximate 



'2A 



h ■i2^N [kptf \eF 

The integral over the spectral function gives 

dr{s^{r,0))' = 2n T dr r(^(r))^ ^ ^J^Mll. 



1/2 



2^ 4 



(7.45) 



(7.46) 



where 



dyye-^y/\J^iy)^^\n2x, x » 1. (7.47) 

The approximate result can be shown with the aid of L' Hospital's rule. Similarly, for a 
screened interaction 



Ui2ir) = C/i2(r) exp(-r/r,), Ui2{r) 



(7.48) 
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Fig. 7.1; Plot of T(x) and 
((x) as a function of x. 
The dashed line shows the 
logarithmic approximation 
to the latter. 



2 4 6 8 10 

X 



we get 



j ArUl^ir.Q) = 2tt dr rUl^i^J r"^ + <P) 



= 2^UUd)d^'^o{'2d/r,) = ^UUd)^ {2d / r ,)rl , (7.49) 

where 

Tc{x) = / dte"^e-\ (7.50) 

J X 

is the incomplete Gamma function and 

T(x) =.T2ro(x). (7.51) 



The two introduced functions C,{x) and T{x) are shown in Fig. 7.1. Collecting things 
we now obtain the estimate 



((5G2i)) :.,.-(^— .— — J — , (7.52) 

where the numerical prefactor is given by 



2-"vW6 - 3.5 X 10"^ (7.53) 
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This result is reported in Ref. [F] which generalizes the result first given in Ref. [E] for 

d ~ r^. 

For typical numbers (such as those in Table 1.1 for the quantum dots shown in 
Fig. 1.5) we estimate the drag fluctuations to be of the order of 0.1 Ohm. As we have 
already noted the zero mean value and the magnitude of the finite fluctuations are not 
changed by breaking of time reversal symmetry, i.e. (3=1 — > 2. The prediction 
of zero mean value and finite fluctations can be considered as a test of the degree of 
ergodicity. 



Chapter 8 



Summary 



In the following we give a brief summary on the mesoscopic drag formalism, the nu- 
merical method, and the applications to quantum wires and quantum dots. Finally, 
concluding remarks are given. 



8.1 Mesoscopic formalism 

Starting from Kubo formalism we have calculated the DC drag conductance G21 which 
relates the frictionally mediated response I2 in subsystem 2 to a driving field Vi in 
subsystem 1. The calculation is based on a second order perturbation expansion in the 
inter-subsystem Coulomb interaction Ui2- The drag conductance. 



2kTsmh^{huj/2kT) 

is expressed in terms of two "triangle" functions = ^{hPiPi}- Whereas often 
Ui2{ri, r2) ~ Ui2{ri — r2) we for disordered systems with broken translation sym- 
metry have Ai(w, r^, r[) 7^ Ai{u!, ri ~ r^). This for instance means that the outcome 
ofEq. ^ can have both signs depending on the given disorder potential. 

Our calculations, which are restricted to Fermi liquids, demonstrate how we with 
the aid of Wick's theorem can formulate the A's in terms of the single-particle prop- 
erties of the two subsystems. The single-particle properties can be expressed by i) the 
wave functions in terms of the matrix elements of the current / and particle-density p, 
ii) by the spectral function Hi) or by scattering states where the current matrix / 
is expressed in terms of the scattering matrix S. Furthermore, we have derived results 
for the low-temperature limit and shown that G21 oc as is also the case for the rate 
for ordinary particle-particle scattering. This formalism provides the framework for 
Refs. [E,FG,H]. 
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8.2 Numerical method 

Using the finite differences approximation we have mapped Eq. ( ^j] ) onto a lattice and 
also formulated how the spectral function can be obtained by use of a simple numerical 
method. The resulting formula for G21 can be written as a trace over matrices and in 
the low temperature regime we get 



e2 -72 . _ 

G21 = -J^ikT^Y {f^2iAi(7i2A2}, (8.2) 

where A is expressed in terms of the spectral function s:/^j, at the Fermi level. Here, 
7 is the hopping element of the tight-binding-like Hamiltonian (band-width is 47). 



Eq. (8.2) is very similar to existing trace formulas for the Landauer conductance Gu. 
The implementation is straight forward and the method provides a very efficient tool for 
investigating drag in different geometries and/or for different disorder configurations. 
In this work the method is applied to ensembles of disordered wires. 



8.3 Quantum wires 

For one-dimensional quantum wires we have shown how in the ballistic regime analyt- 
ical progress can be made. For long wires, kpL ^ 1 we find that 

Gai oc t/f2(2fcF), e — ^00, (8.3) 

where Urziq) is the Fourier transformed interaction and i the mean free path. For 
quasi-ballistic wires the corrections to the mean-value are small so that 



(G21) oc C/f2(2fcF) + Oil/kpe). (8.4) 

For the fluctuations we apply perturbation theory and to leading order in the disorder 
strength we find that 

{{SG2if) cx {TZyui2{2kF)UU0), (8.5) 

where (TV) ~ L/£ ^ 1 is the mean reflection probability of the wires. For long-range 
interaction the relative magnitude 



(G20 ^^"^^iwk^y ^'-'^ 

can be of the order of unity because C/i2(0) ^ Ui2{2kF). The large fluctuations 
can be understood from a combination of back-scattering in the disorder channel and 
forward-scattering in the Coulomb channel — breaking the translational invariance by 
disorder allows for transfered momenta different from 2fci? and this gives rise to large 
drag fluctuations [E]. 
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This we have also studied numerically and we found excellent agreement with the 
perturbative results [E] . In addition to the statistical properties the numerical studies 
have also given insight about the full distribution of G21 for different natures of the 
disorder [F, G, H]. If the disorder potentials of the two wires are mutually fully cor- 
related we diagrammatically find that mean value and fluctuations are enhanced com- 
pared to the un-correlated case, for the latter by a factor %/2- These predictions are also 
confirmed numerically [F, G, H] . 



8.4 Quantum dots 

For coupled chaotic quantum dots we have used random matrix theory to calculate 
the statistical properties. For the mean value we find that there is no drag on average, 
(G21) = 0. The physical explanation is that in the chaotic regime the induced current 
has no preferred direction and electrons will with equal probabilities escape the dot 
through either of its two leads. This implies finite fluctuations which is also what 
random matrix theory predicts. The estimate that we get for the fluctuations is of the 
form [E, F] 



.2\i/2 krVA 
{i^G2i) ) oc , (8.7) 

where the screening length is of the order of the separation d, A is the interaction 
area, and N the number of modes. Using realistic numbers we estimate that the drag 
fluctuations can be of the order of 0.1 Ohm. 



8.5 Concluding remarks 

Rojo ends his 1998 review on electron-drag [ ]66[ | by the sentence: "The magnitude of 
die chaUenge seems to equal that of the progress made so far in the field". Since then 
the study of the mesoscopic regime has been initiated [TI] , E] and though there has 
been substantial theoretical progress I belive the major challenge is on the experimental 
side. However, the results reported in this thesis and attached publications [E, F, G, H] 
and by others [ p9| , [7T| ] illustrate the potential reward of overcoming the experimental 
challenges: Although the mesoscopic fluctuations of the Coulomb drag conductance 
are non-universal they are an extreme example of mesoscopic fluctuations — they can 
be of the order of, or even exceed, the mean value! Probably other mesoscopic phe- 
nomena are to be discovered and I hope that this work may serve as an inspiration for 
experimentalists to conduct Coulomb drag measurements in the mesoscopic regime. 
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